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Abstract 



■ Quantum computers hold the promise to solve certain computational task much 

^ i more efficiently than classical computers. We review the recent experimental ad- 

vancements towards a quantum computer with trapped ions. In particular, various 
implementations of qubits, quantum gates and some key experiments are discussed. 
Furthermore, we review some implementations of quantum algorithms such as a 
OO ' deterministic teleportation of quantum information and an error correction scheme. 
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1 Introduction 



The aim of tliis article is to review the recent development of ion trap quantum 
computing. The field evolved rapidly in the recent decade. Thus the many 
facets of experimental ion trap quantum computing and its techniques cannot 
be covered all. Instead, we want to present here a coherent picture of the most 
important experimental issues and refer the reader to the original publications 
for the details. We also describe some of the milestones achieved so far in ion 
trap quantum computing, like teleportation of quantum states and quantum 
error correction. Still much of the work especially towards shuttling ions with 
segmented traps is only touched upon. 



A quantum computer uses the principles of quantum mechanics to solve cer- 
tain mathematical problems faster than normal computers. Such a quantum 
computer processes quantum information whose most basic unit is called a 
quantum bit (qubit). Already a small quantum computer, consisting of forty 
qubitJ^. could solve quantum mechanical problems that are intractable with 
current computers. In particular, the study of quantum mechanical many body 
syste ms would benefit considerably from such a device (IFeynmanl . Il982l : iLloydl . 



19961 ). In 1989, David Deutsch discovered a mathematical problem which can 



be solved faste r with quantum mechanical means than with classical ones 
( jPeutschl . Il989l ) . But it was a few years later when the rapid development of 
quantum computation set in, marked by Peter Shor's discovery of a quantum 
algorithm with which large numbers ca n be factored much faster than with 
today's classical algorithms (IShorl . 119941 ). 



Shortly afterwards, Ignacio Cirac and Peter Zoller found a physical systern 
on w hich such quantum algorithms could be implemented (ICirac and Zoller . 



19951 ): single trapped ions were supposed to carry the quantum information. 



which is manipulated and read out with focused laser beams. Already within 
a year's time, David Wineland's group at National Institute of Standards and 
Technology de r nonstra ted the heart of such an ion-trap quantum computer 
( iMonroe et al.l . Il995al ). a controlled bit flip on a single ion. Without exagger- 
ation one can say that those two publications mark the birth of experimental 
quantum computation. Rapidly also other implementations were considered. 
In particular, liquid-sta te nuclear magnetic reson a nce was used to demonstrate 
a quantum algorithm (IGershenfeld and Chuangl . 119971 : IChuang et al.l . Il998b : 



Jones and Moscal . Il998f) and later even Shor's a lgorithm was implemented on a 



seven qubit register (jVandersypen et al.l . 120011 ). On the theoretical side, many 
new interesting implications were found. For quantum computation, the most 



^ To describe the (arbitrary) state of a forty qubit system, 2^^^ complex numbers 
are necessary. Already this requirement exceeds the capacities of current super com- 
puters. 
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relevant impli cation was t he discovery of quan tum error cor rection protocols 
by Peter Shor (ISlioii . ll995l ) and Andrew Steane (ISteand . 119961 ) . These protocols 
allow for the implementation of arbitrary long quantum algorithms without 
perfect control. 

The basic unit of quantum information (a qubit) can be implemented with 
a two-level system such as the electron's spin in a magnetic field or using 
two levels of an atom. A simple quantum computation initializes the qubits, 
manipulates them and finally reads out the final states. Any physical imple- 
mentation of quantum computation must be able to perform these tasks. Thus 
the physical system must satisfy th e requiremen t s laid out in Sec. l2.1l to qualify 
as a universal quantum computer (jDiVincenzd . l200ll ) . 



(1) It must be a scalable physical system with well characterized qubits. 

(2) It must be possible to initialize the qubits. 

(3) The qubits must have a coherence time much longer than the operation 
time. 

(4) There has to be a universal set of quantum gates. In the most simple case 
one considers single-qubit and entangling two-qubit gates. 

(5) A qubit-specific measurement must be attainable. 

Furthermore, to build up a quantum network, one requires: 

(6) The ability to interconvert stationary and flying qubits. 

(7) The ability to faithfully transmit flying qubits between specified locations. 

In principle, these requirements can be fulfilled with a number of physi- 
cal a pproaches, like nuclear magnetic resonance (IGershenfeld and Chuand . 



19971 ). cavity quantu m electrodynan i ics (c avity-QED) (IRaimond et al.l . l200ll ) 



Josephson junctions (IMakhlin et al.U200l[) . a combination of ci rcuit-QED and 



Jose phson junctions (IWallraff et al.l. |2004 iMaier et al.l . 120071 ) and quantum 



tics have also been proposed (IKnill et al.l. l2001b ) and actually were used fo r 



dots (iLoss and DiVincenzol . ll998UPetta et al.l.l2005h . Methods using linear op 



quantum information processing (IWalther et al.l . l2005l : iLanyon et al.l . l2008l ). 
however, here the emphasis is naturally more on the transmission of quantum 
information rather than processing it. 



Using nuclear magnetic resonance, quite a number of impressive demonstra- 
tion experiments have been performed. Unfortunately, the state of the quan- 
tum re g ister (molecules) can only be poorly initialized (except in special cases 



(jjonesl . |2000|)), making t he scaliri g properties o f NMR quantum computa- 
tion not very promising flWarrenl . Il997i : Ijonesl . I2OOOI : iLinden and Popescul 



2OOII ). At the moment Josephson junctions seem to be very a ppealing, espe 



daily when combined w ith superco nducting strip line cavit ies (IWallraff et al 



2004 



Majer et al.l . 120071 ). Recently. ISteffen et al.l (120061 ) and lPlantenberg et al. 



(I2OO7I ) demonstrated Bell-states and a controUed-NOT, respectively, with pairs 
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of Josephson junction qubits. 



2 Ion trap quantum computers 



Of the many approaches that have been proposed for construc ting a quantum 



computer, trapped ions are currently one of the most advanced (lArmy Research Office (USA) 



20051 ) ■ This is also reflected in the fact that around the year 2000 only half 



a dozen groups pursued experimental ion-trap quantum computing while in 
2008 there are more than 25 groups working in the fleld. 

Already long before the idea of quantum computation was picked up by exper- 
imentalists, four out of the flve core criteria required by DiVincenzo for a quan- 
tum cor aputer were denaonstra ted with trapped ions in the laboratory: initial- 
ization (Wineland et al.. 1980) and read-out of the internal electronic states 



of tra pped ions (INagourney et al.l . Il986l: ISauter et al.l. 1 1986 



19861 ) , extremely long coher ence times (IBollinger et al.l. Il991l) and laser cooled 



ion crystals with many ions (jPiedrich et aL 



1987 



Bergquist et al. 



Wineland et al.l . 119871 : iRaizen et al. 



Il992a .b.) serving as a qubit register. Then lCirac and ZoUeri (il995l ) realized that 
a quantum computation can be carried out by couphng the ions via a collective 
motional degree of freedom. Thus, a route to implement the missing condi- 
tional evolution of physically separated qubits (a two-qubit interaction) was 
introduced. In addition they proved that the size of the resources necessary 
to cont rol trapped ions does n ot increase exponentially with the number of 
qubits ( Cirac and ZoUer . 19951 ). 



Soon after the proposal by Cirac and Zoller in 1995, the NIST ion stor- 
age group around David Wineland implemented the key id ea of the pro- 
posal — a conditioned phase shift — with a single Be^-ion (IMonroe et al. 



199531). Furthermore they demonstra t ed a few other two-qubit gate candidates 
flSackett et aP . I2OOOI: IPeMarco et all I2OO2I : Iheibfried etaP . l2003bh , entangled 



up to four io ns (ISackett et al 



free subspace (iKielpinski et al. 



flLeibfried et al.l . l2002h . 



2000h . demonstrated a so-called decoherence 



2OOII ) and simulated a nonlinear beam-splitter 



In our group in Inns bruck, the Deutsch -Josza algorithm was demonstrated 
with a single Ca"'"-ion (I Guide et al.l . 120031 ) , follow ed by the flrst implernentatiq n 



of a set of universal gates on a two-ion string (jSchmidt-Kaler et al.l . l2003d ) 



In addition, the creation of various entangled states and the par t ial rea d-out 
of an entangled quantum register was demonstrated flRoos et al.l . l2004bl ). 



Further milestones in ion-trap qua ntum computing were experiments on 



tum teleportation by both groups (iBarrett et al.l. |2004J: iRiebe et al.l. |200 



error correction protocol by the NIST-group (IChiaverini et al.l . |2004J ). entan- 



uan- 
an 
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glement of six and eight particles (ILeibfried et all l2005l : iHaffner et all l2005a ) 
by bo th groups, respectively and entanglement purification ( iReichle et al. 



2006bl ). Recently,ions in separate traps have been also entangled using ion- 



photon entanglement by the Ann- Arbor group (IBlinov et al.l . l2004l : iMaunz et al. 



20071 : iMoehring et~aD . l2007l : iMatsukevich et al.l . [2008h 



Currently, miniaturization and integration of segmented ion traps is rapidly 
progressing. Already for some time, the NIST group has been successfully us- 
ing microfabricated segmented traps to relieve the difficulties of single ion ad- 
dressing. Inspired by this success, in the mean time virtually all ion trap groups 
started to develop segmented trap technologies. Furthermore, U.S. funding 
bodies initiated contacts between the various groups to further ion trap re- 
lated technologies and established contacts to microfabrication laboratories, 
such as Lucent Technologies and Sandia National Laboratories. In Europe the 
'Specific Targeted Research Project' MICROTRAPS with six participating 
groups has been formed to further microfabricated ion trap technologies. 



2.1 Principles of ion-trap quantum computers 



A excellent overview and detail ed account of the fundamental issues of ion trap 
quant um computing i s given bvlWineland et al.l (119981 ) and by lSasura and Buzek 
( I2OO2I ). Furthermore. ILeibfried et al.l ( 2003al ) review the progress towards the 
manipulation and control of single ions. Very recently, the genera t ion a nd 
applications of entangled ions were reviewed by iBlatt and Winelandl fl2008h . 



We start here by summarizing how i on-trap quantum computers fulfill the 
DiVincenzo criteria mentioned above (iDiVincenzol . 1200 ll ): 



(2) 



(3) 



A scalable physical system with well characterized qubits: long- 
lived internal levels of the ions serve as the qubits (see Sec. 12.31) . The 
qubit register is formed by strings of ions in a (linear) trap. While this 
approach is in principl e scalable, it is desira b le to distribute t he io ns 
among multiple traps ( Wineland et al. . 19981 : Kielpinski et al. . 2002)). 
Thus complications, due to the increased mass and more the complicated 
mode structure of large ion strings can be circumvented. First steps in 
this direction will be briefiy discussed in Sec. [HI 

The ability to initialize the state of the qubits: this is most easily 
achieved by optical pumping to a well-defined electronic state. Fideli- 
ties of 0.99 are typical, in Sec. 12.41 also methods to achieve much higher 
fidelities will be discussed. 

A coherence time much longer than the operational time: in cur- 
rent quantum computing experiments, typically coherence times of a few 
milliseconds are achieved which are about one to two orders of magni- 
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tude longer than the time scale for quantum operations (see Sec. [3] for 
details). The coherence time is often limited by magnetic field fluctua- 
tions, on the other hand coherence times of more than 10 s have been 
demonstrat ed with Raman-tran sitions between magnetic-field insensitive 
transitions (ILanger et al.l . l2005l ). Simil ar observations were made in ^^Ca ^ 



(4) 



by the Oxford and Im isbruck groups (iLucas et al.l . 120071 : iBenhelm et al. 



2008al ). Furthermore, iBoUinger et al.l (Il99ll ) demonstrated a coherence 



time of more than 10 minutes using a microwave drive instead of a Raman- 
laser setup. However, so far only one experiment demonstrated the com- 
bination of high-fidelity quantu m gates and magnetic field insensitive 
transitions (IHaljan et al.l . l2005bl ). 
A universal set of quantum gates: 



(a) 



Single qubit gates are implemented by driving Rabi oscillations be- 
tween the two qubit levels with resonant laser pulses (see Sec. 12. 5p . 
The gates can be represented as rotations of the Bloch sphere where 
the axis of rotation can be selected by changing the phase of the ex- 
citing laser field (single photon transition) or the phase difference of 
the two Raman beams (see Sec. 12.51) . A single-qubit phase-gate can 



be produced directly by an off-resonant laser via an AC-Stark shift 
(see Sec. 12.51) . 

(b) For two-qubit gates usually the long range interaction due to the 
Coulomb force is employed. In the original proposal by Cirac and 
ZoUer, the quantum information of one ion is swapped to the commo n 
motional degree of freedom of the ion string (jCirac and ZoUen . Il995l ). 
Then an operation conditioned on the motional state can be carried 
out on a second ion before the quantum information is swapped back 
from the motion to the first ion. Section 12.61 details this idea as well 
as other methods to implement multi-qubit gates. Implementations 
of some of the two-qubit gate recipes will be discussed in Sees. 14.11 
and 1121 

(5) A qubit-specific measurement: one of the qubit levels is excited on 
a strong transition to a higher lying auxiliary short-lived level while the 
other qubit level remains untouched. Thus, fluorescence from the decay is 
detected only if the qubit is projected to the qubit level which is coupled 
to the auxiliary transition (see Sec. 



Additionally DiVincenzo requires: 



(6) 



(7) 



The ability to interconvert stationary and flying qubits: ions can 
be stored in a high-finesse cav ity. Thus, the ions ' internal state can be 



mapped onto a photonic state (j Cirac et al.l . 119971 ). 



The ability to faithfully transmit flying qubits between specified 
locations: a photon can be transmitted through a fiber and at the target 
location coupled via another high-finesse cavity to the target ion. 
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|0,g> 

Fig. 1. Energy level scheme of a single trapped ion with a ground {\g)) and an excited 
(|e)) level in a harmonic trap (oscillator states are labeled |0), |1), |2), • • • ). denotes 
the carrier Rabi frequency. The Rabi frequency on the blue sideband transition 
|0,e) ll,^) transition is reduced by the Lamb-Dicke factor r] as compared to the 
carrier transition (see Eq. [5]). The symbols iOquUt and tot denote the qubit and the 
trap frequency, respectively. 

There exist further possibilities to connect two distant ion trap quantum com- 
puters. One of them will be briefly discussed in Sec. 15.51 

For practical applications, the DiVincenzo criteria have not only to be fulfilled, 
but also the fidelity and speed of the implementations have to be considered 
(see Sec. [8]). Furthermore, it is highly desirable to implement all operations as 
parallel as possible. 



2.2 The basic Hamiltonian 



We will now briefly discuss the Hamiltonian of two-level systems interacting 
with a quantized harmo nic oscillator via l aser l ight. For more detaile d discus - 
sions, we refer to Refs. (IWineland et al.l . Il998l ) and (ILeibfried et al.l . l2003al ). 
The basic level scheme of such a system is displayed in Fig. [H We start out by 
writing the Hamiltonian for a trapped single ion interacting with near resonant 
laser light, taking into account only two levels of the ion and one vibrational 
mode (taken to be oriented along the z direction): 



exp [17] 



ae 



+ h.c. . 



Here, a± is either the atomic raising or the atomic lowering operator, while 
a^ and a denote the creation and annihilation operator for a motional quan- 
tum, respectively. Q characterizes the strength of the laser field in terms of 
the so-called Rabi frequency, ip denotes the phase of the field with respect to 
the atomic polarization and A is the laser-atom detuning. Ut denotes the trap 
frequency, t] = k^zo is the Lamb-Dicke parameter with kz being the projection 
of the laser field's wavevector along the z direction and Zq = Jh/ {2miUt) is 
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the spatial extension of the ion's ground state wave function in the harmonic 
oscillator (m is here the ion's mass). We mention also that the rotating wave 
approximation has been applied which assumes that both the laser detuning 
and Rabi frequency are much smaller than optical frequencies. A similar treat- 
ment can be carried out for qubits based on Raman-transitions by eliminating 
the virtual level through which the two qubits are coupled. Please note that 
in our definition the Rabi frequency measures the frequency with which the 



popul ation is exchanged iri contra st to the definition used by IWineland et al 



fll998f ) and lLeibfried et all fl2003af ) 



Using the Lamb-Dicke approximation [rjJ {{a + a^Y) <^ 1) whic h is almost al- 



ways valid for cold tightly bou nd ion strings, we can rewrite Eq. [2] (iLeibfried et al 
2nn8al : I.Tonathan et"aD . boonh : 



H = m[a+e-'^^'-^^ + a.e'^^'-^^ (2) 

Three cases of the laser detuning A are of particular interest (see Fig.[T]): A = 
and A = icut. This becomes apparent if a second rotating wave approximation 
is carried out where it is assumed that only one transition is relevant at a time. 
Discarding time dependent terms, we thus arrive at 

(1) the Hamiltonian describing the carrier transition (A = 0) 

H'^''' = hQ{a+e'^ + a.e"^^) . (3) 

Here only the electronic states \g) and |e) of the ion are changed. 

(2) the Hamiltonian describing the blue sideband transition (A = Ut) 

= iWL7]{a+a}e''^ - a^ae-'"^) . (4) 

Simultaneously to exciting the electronic state of the ion, in this case a 
motional quantum (a phonon) is created. Within this two-level system, 
Rabi flopping with the Rabi frequency 



VLn,n+i = Vn + lr] (5) 

occurs, where n describes the number of motional quanta (phonons). For 
convenience, we define the blue sideband Rabi frequency f2+ = f2o,i which 
describes the flopping frequency between the |5',0) and the |e, 1) state. 
(3) the Hamiltonian describing the red sideband transition (A = —Ut) 

H- = ihQr]{a_a^e-'^ + a+ae''^) . (6) 

Simultaneously to exciting the electronic state, here a phonon is destroyed 
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and Rabi flopping with the Rabi frequency 



^n,n-l = y/uT] Q (7) 

takes place. 

Naturally, one is interested in driving the sideband transitions as fast as pos- 
sible to speed up the quantum operations based on the sideband transitions. 
However, especially for small Lamb-Dicke factors, the second rotating wave 
approximation performed to obtain Eqs. [3ll6] is then not satisfied for Rabi fre- 
quencies Q comparable to the trap frequency Ut : as can be seen from Eq. [5] for 
ions in the motional ground state, the carrier transition is stronger by a factor 
of as compared to the sideband transitions. Therefore, driving the weak 
sideband trans itions with strong la ser fields. Stark shifts and off-resonant ex- 



citations arise (ISteane et al.l . l2000l ). For example, to achieve a side-band Rabi 
frequency fl^ of a fraction / of the trap frequency uJt, we need a laser field of 
strength Q = ^l+frj = fut/r]. The Stark shift AE of the qubit transition due 
to the presence of the carrier transition is given by 

AE n"^ 1 fu;t^+ f 



h 2A 2cJt 2r]^ 



The phase evolution AEt/h due to the AC-Stark shift becomes comparable 
to the desired Rabi flopping with frequency Qj^ already for f ^ rf and thus 
typically for / ~ 0.01 [t] ~ 0.1). Similarly, off-resonant excitations on the 
carrier transition might spoil the fidelity in this regime. In pa rticular, Rabi 



oscill ations on the carrier transition occur with amplitude A ( ISteane et al. 
2000h 



^ 09 , . .2 ~ , ,2 „2 ' 



where we assumed Q Ut, which is justified by the conclusions from Eq. [HI 
Both the AC-Stark shift and the off-resonant excitations, can be at least par- 
tially canceled with methods described in Sec. 13.1.31 However, it still remains 
very difficult to drive Rabi flops on sidebands much faster than rjUf Thus, it 
is hard to implement a two-qubit gate based on sideband transitions within 
one trap period. 



2.3 Choice of qubit ions 



A good qubit candidate must meet certain criteria. In the case of trapped 
ions it is usually sufficient to concentrate on long coherence times on the 
qubit transition compared with the manipulation times and on the technical 
feasibility of the required lasers. The scale to which the coherence time has to 
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Fig. 2. Currently achieved time scales in ion trap quantum computing. All operations 
can be implemented faster than the relevant decoherence mechanisms. 




Fig. 3. Generic level schemes of atoms for optical qubits (left) and radio-frequency 
qubits (right). In addition to the two qubit levels {|0), |1)} usually a third rapidly 
decaying level is used for laser cooling and state read-out. While the optical qubit is 
typically manipulated on a quadrupole transition, the radio-frequency qubit levels 
are connected with Raman-transitions. 



be compared with is given by the gate operation times of typically 0.1 

(c.f. Fig ED. 



'500/is 



Currently two different scliemes are used to store quantum information in 
trapped ions: In the first scheme, superpositions between the electronic ground 
state and an excited metastable electronic state provide the two-level system 
(optical qubit) (see Figs. [3^ andH]). The excited D.^/^-Ieve l in Ca"*" — similarly 



to Sr"*" and Ba " ^ — has a life time of more than one second (IBarton et al.l . 12000 



Kreuter et al.l . l2004l . 120051 ). In a second scheme, even larger coherence times 
can be achieved with superpositions encoded in the electronic ground-state of 
the ions (radio- frequency qubits) (see Fig. [3]). Here, either the Zeeman or the 
hyperfine structure can be used. The lifetimes of these states are estimated to 
be much larger than any currently experimentally relevant timescales. 
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There are different aspects to be considered for optical and and radio-frequency 
qubits: to fully exploit the potential of the optical qubit provided by the 5*1/2 
and D^/2-\eve\ in Ca"*", a laser line width of less than 200 mHz is required. This 
fractional stability of about 10~^^ is technically quite demanding, however, it 



has alread y been implemented for metrology purposes flYoung et all ll999al Jbl: 



Hall et al.l . |2006| ). The transition frequencies for radio- frequency qubits are 



usually below 10 GHz. Thus it is much easier to get the required stable phase 
reference. 



Radio-frequency qubits can be manipulated either directly by microwaves 
( iMintert and Wunderlichl. 120011) or — as it is usually the case — by a Raman- 
process ( iMonroe et al.l . Il995al ). Here, a single atom absorbs a photon from one 
laser field and emits a second one into a second laser field. The frequency dif- 
ference of the two laser fields provides the necessary energy for the population 
transfer and thus only this frequency difference is important. The two laser 
fields are either derived from the same laser or a phase locking between two 
lasers can be implemented without resorting to optical frequency standards. 
Using optical fields instead of microwaves has the advantage that a much 
better spatial concentration of the power can be achieved and much higher 
Rabi-frequencies can be attained as with microwaves. Furthermore, by driving 
radio-frequency qubits optically with two anti-parallel beams, two photons can 
transfer their recoil onto the atom. Thus, the coupling to the ion motion is 
increased as compared to single photon transitions (larger Lamb-Dicke factors 
7]) and as a consequence the speed of two-qubit gates can be higher. Choosing 
co-propagating beams, the coupling to the motion can be inhibited efficiently 
which has the benefit of suppressing the sensitivity to the ion motion. 



For optical as well as for rf-qubits, the coherence times are often limited by 
fluctuations of the magnetic fields. The reason is that usually the qubit basis 
states employed have different magnetic moments such that they experience an 
additional phase evolution due to the (fluctuating) magnetic field. Strategies 
to avoid these decoherence sources will be discussed in Sec. I3.1.2[ 



We note here that a good qubit must not necessarily combine a large co- 
herence time to manipulation time ratio with high fidelity initialization and 
read-out capabilities: initialization and read-out can be implemented with an 
additional auxiliary ion of a differe nt ion species. This id ea has already been 
demonstrated by the NIST group (jSchmidt et al.l . |2005| ): In this experiment 
aimed at implementing a frequency standard, the state of an individual A1+ 
ion was detected via a Be"*" ion. However, for practical purposes it is desirable 
to initialize and read out the qubit ion directly. 



Finally, a good qubit candidate must have all the relevant transitions in an 
accessible frequency regime. Generally, laser sources, fibers and detectors for 
short wavelengths are more expensive, less efficient and often more fault-prone 
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than for longer wavelengths. However, ions tend to have short wavelength 
transitions. Thus there are only a very limited number of ions which have 
strong transitions in the visible frequency range. In the quantum computing 
context, calcium, strontium and ytterbium ions appear to be attractive due to 
their relative large wavelength transitions. However, beryllium and magnesium 
have a relatively small atomic mass which leads to large Lamb-Dicke factors 
(?7 ~ 0.3 in some of the NIST experiments) and eases coupling of the electronic 
and motional degrees of freedom and thus makes them attractive in spite of 
difficulties with laser radiation generation, manipulation and fiber optics. 



2.4 Initialization and read-out 



Prior to the implementation of a quantum algorithm, the qubits must be 
initialized in a well-defined state. For atoms, in general, this can be m ost con- 
venien tly achieved by optical p umping. This idea w as intr oduced bv iKastlei 
( 



1977h and flWineland et al.LllQSd ). 



1950h and first implemented by Brossel et al. ( 1952 ) and by Hawkins and Dicke 



1951 . Fo r reviews on optical p umping see Refs. flHappeil . Il972f ). dWebeil . 



The general idea of optical pumping is that an atom is driven until it decays 
into a state where the drive does not act any longer. Usually circularly polar- 
ized light is used to pump the atom into one of the extreme Zeeman-levels. 
Typically, the target state is occupied in less than 1 /xs with a probability 
larger than 0.99. The initialization fidelity is usually limited by the quality of 
the polarization of the driving laser along the preferred axis of the qubit. In 
most cases, this axis is given by the direction of the magnetic field. 



For fault-tolerant quantum computing, however, it appears that initialization 
fidelities exceeding 0.9999 are desirable. It is not clear whether such high 
preparation fidelities can be achieved with the methods described above. Fur- 
thermore, there exist situations where it is not possible to achieve a pure 
circular polarization along the quantization axis. The latter might happen for 
instance, because there is no optical access along the direction of the magnetic 
field to send a laser in. In these frequency rather than a polarization 

selection c an be used for opt ical pumping if the ions offer a spectrally narrow 



transition (IRoos et al.l . l2006l ). 



We illustrate this procedure here for the Ca"*"-ion: with a narrow band laser 



the ion is excited on the S1/2 



+1/2) 



D5/2(w^j = —3/2) transition 



(see Fig. H]) and simultaneously the D level is coupled to the P3/2 level with a 

TUj = +1/2) 

3/2 [mj = {-3/2,-1/2}) 



broad band laser at 854 nm. In this way the population of the S1/2 1 
state is effectively coupled to the rapidly decaying P 



levels while the S1/2 



+ 1/2) level is not touched. Fig. [S] shows the deple- 
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T ~ 1.2 s 



m.= + l/2 



m:=-l/2 



Fig. 4. Level scheme of ^''Ca"*", with Zeeman substructure and required laser wave- 
lengths for manipulation of the calcium ions. 



tion of the initially fully populated S 



1/2 



-1/2) level with time. Pumping 



time constants of about 10 /is have been demonstrated in Innsbruck. 

In principle, the attainable pumping efficiency is only limited by the strongest 
Rabi frequency and pumping time. Assuming a typical transition splitting of 
the various Zeeman transitions of a few MHz, Rabi frequencies ^Rabi lower 
than 27T x 10 kHz are required to keep the off-resonant excitations Pofr ~ 
f2^/A^ (see Eq. [9]) of the unwanted transitions below 10~^. From that, we 
estimate that within 1 ms an initialization fidelity of 0.9999 can be reached. In 
experiments in Innsbruck, prepar ation fidelities exceeding 0.999 have already 
been observed (IRoos et al.l . l2006l ) . 



At the end of a quantum algorithm the quantum register needs to measured. 
Radiation couplin g to only one of the qubit levels can be used here. This idea 
was introduced by iDehmeltl (1l975l ) and is often termed electron shelving. Ex- 
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Fig. 5. Optical pumping of the Ca"*" S1/2 {mj = +1/2) level using frequency 
selection rather than polarization selection. 



perimentally it was first implem e nted by iNagourney et al.l (ll986l ). ISauter et al. 
fll986h . and bv iBergquist et all fll986h . 



Referring to the level scheme of ^''Ca"'" in Fig. HI the ion does not fluoresce 
under irradiation of light on the S1/2 — > P1/2 and 1^3/2 —>■ P1/2 transitions if its 
valence electron is in the D5/2-state. If the electron, however, is in either the 
S1/2-, Pi/2- or -D3/2-state, the ion will scatter approximately 10^ — 10^ pho- 



tons/s. A lens system collects typically 10~^-10^^ of this fluorescence light 
such that with a photomultiplier tube (quantum efficiency ~ 30%) about 
30 photons/ms can be detected. As typical background count rates of photo- 
multiplier tubes are usually well below 1 count/ms, one expects exceedingly 
high state detection fidelities when collecting fluorescence for detection times 
of more than one millisecond. However, this reasoning holds only if the ion 
has a negligible probability to decay either from the P1/2 to the 1^5/2 or from 
the D5/2 to the 5*1/2 or the D3/2 level during the detection time. In practice, 
these two constraints lead to an optimal detection time where the error due to 
the Poissonian spread of the number of scattered photons is balanced with th e 
relaxation time scales (Roos, 2000l : Acton et al. . 2006 : Myerson et al. . 2008). 



Typical detection times are about one millisecond with detecion fidelities ex- 
ceeding 0.99. The fidelity can be further increased by anaysing the arrival 
times of the photons with a maximum likelihood method and thus identifying 



some events when the 1)5/2 level decayed during the detection (IMyerson et al. 



20081) 



Other qubits can be detected very similarly. If the energy separation of the 
two qubits, however, is not large enough to allow for selection via the laser 
frequency, the polarisation of the laser field can be used. For instance for 
^Be"*", circular polarization ensures that predominantly only one of the two 
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qubit states scatters photons ( Sackett et al. . 20001 : Langer . 2006 ). 



The detection efficiency can be further increased when the quantum state is 
mapped onto auxihary qubits prior to its measurement: for this, one prepares 
first an additional qubit (the ancilla qubit) in state |0). A controUed-NOT op- 
eration maps then the qubit a|0)+/?|l) onto the combined state a|00)+/3|ll) of 
the two qubits. Measuring both qubits yields an improved fidelity as compared 
to measuring a single qubit if the fidelity of the ancilla qubit state prepara- 
tion and the controlled-NOT operation are high enough. ISchaetz et al. 1 (l2005h 
demonstrated this procedure, albeit preparing the ancilla ion in (|0) + \ \))/\/2 
and replacing the contolled-NOT with a controlled phase gate (see Sec. 12.6.41) 
followed by a single-qubit rotation on the ancilla. 



Hume et al.l (120071 ) present another variant to improve the state detection 
fidelity by measuring the qubit repeatedly. In this experiment, the qubit is 
encoded in superpositions of the ground state and the excited qubit level of 
^^Al"*". This qubit state is fi rst transferred to a ^Be"*" ion by a series of laser 
pulses (jSchmidt et al.l . l2005l ) : a ffist pulse on the red sideband on an auxiliary 
transition of the Al"*" ion inserts one motional quantum only if the A1+ ion 
is in the ground state. This motional excitation is transferred to the ^Be^ 
ion with a red sideband pulse on the ^Be"*" ion which is then detected via the 
usual state dependent fluorescence method. Most importantly, the electronic 
auxiliary state of the Al"*" ion decays back to the ground state on a timescale 
of about 300 /xs and not to the excited qubit state. Thus the information 
in which state the qubit is projected is still available in the A1+ ion and this 
procedure can be repeated. Based on the photon counts from the ^Be+ ion and 
previous detection results, another d etection step is ca rried out. Repeating 
this procedure on average 6.5 times, iHume et al.l (120071 ) achieve a detection 
efficiency of 0.9994. The accuracy of this method is only limited by processes 
which couple the excited qubit state of the A1+ ion (~ 21 s) to one of the 
electronic levels which does not decay to the ground state. 



Very recently, iMyerson et al.l (120081 ) took adavantage of the arrival times of 
the photons and combined it with the adaptive scheme just discussed. Thus 
they reached within an average detection time of 145 /is efficiencies of about 
0.9999 of a qubit encoded in the S — D manifold of a single '^'^Ca"''. 



2.5 Single-qubit gates 



It can be shown that all quantum algorithms can be broken down into a 
sequence of single-qubit operations plus a specific two-qubit operati on, e.g. a 
cond it ional phase gate, a control led-NOT gate or a VSWAP gate (IDeutschl . 
19891 : iNielsen and Chuangl . I2OOOI ). In ion traps, single-qubit operations are 
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usually easy to implement and thus it is reasonable to use this approach to 
attain universal quantum computing. 



Rabi oscillations between the two qubit levels (see Fig. [6]) implement such 
single-qubit operations. Mathematically, we can describe the effect of resonant 
radiation inducing suc h a coupling by a rotation if ^jO, ip) acting on the state 
vector a|0) + (c.f. iNielsen and Chuand fl2000h ): 



R 



Ci 



exp [iQ 1 2 (^e*'^cr+ + e *'^cr_ ^ ^ = / cos 0/2 + z (cXj; cos Lp — ay sin Lp) sin 6/2 

(10) 

ie"*^ sin 0/2 cos 0/2 



where a+ 



1 




and (T- 





1 



are the respective atomic raising and 



lowering operators, 



1 



and (T,, 



-i 



are the Pauli-spin matri- 



1 oy \i Q 

ces. The angles and ip specify the rotation. 



Often single-qubit operations are visualized by use of the so-called Bloch 
sphere. We identify the north pole of the Bloch sphere with logical |0) (the 
energetically higher state) and the south pole with |1) (see Fig. [7j). In the 
Bloch picture, the angle ip specifies the axis of rotation in the equatorial plane 
and the rotation angle (pulse area), and thus any linear combination of cr^. 
and ay operations can be implemented with laser pulses. 

Rotations around the z axis can be decomposed into rotations around the 
X and the y axis. Alternatively, all following rotations on that qubit can be 
shifted by — A<^ to achieve effectively a rotation about the z axis by /S.(p. 
Finally, a far detuned laser beam can shift the relative energy AE of the 
eigenstates due to an AC-Stark effect by Ai? = r2^/2A (c.f. EqlH]). Thus, after 
a time t = hAip/AE the desired phase shift is acquired. 

The relevant control parameters in the ion trap experiments are the pulse area 
= Qt (Rabi frequency Q, pulse length r) and the phase of the laser field 
ip>. These control parameters can be conveniently co ntrolled with an aco usto- 
optical modulator in double-pass configuration (c.f. iDonley et al.l (120051 )) by 
changing the amplitude and phase of the RF field driving the acousto-optical 
modulator. 



We will give here a quick interpretation of the phase (p: After optical pumping 
the experiments start with the ion qubits in an energy eigenstate. Thus the 
electric field of the resonant excitation laser builds up a dipole (or quadrupole) 
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Fig. 7. Rotation around the y axis visualized on the Bloch-sphere. 

moment oscillating in phase with the field at the laser frequency. In this way 
the first laser pulse (whose length is not a multiple of vr) sets the phase refer- 
ence for all subsequent operations on that ion. Thus, it becomes very intuitive 
to see e.g. that shifting the phase of a second excitation field by vr inverses the 
evolution. 



In ion traps single- qubit manipulatio ns are routinely carried out with fideli- 
ties exceeding 0.99 (iKnill et al.l . l2008l ). Single-qubit-gate fidelities are usually 



limited by laser intensity fluctuations and in the case of single photon tran- 
sitions by the finite temperature of the ion crystal and in the case of Raman 
transitions also by spontaneous emission from the intermediate level. In the 
Lamb-Dicke limit (i.e. the extension of the ground-state wave function is much 
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smaller than the projection of wavelength of the light onto the motion), the 



effect ive Rabi frequency f2efr of the transition is given by (jWineland et al. 



19981) 



n,s^nY[{i-n,7]^), (11) 

i 

where i labels all motional modes, rii is their vibrational quantum number and 
r]i denotes the corresponding Lamb-Dicke factor. Here, in contrast to Eq. [21 
the second order in rj was taken into account. Q is the Rabi frequency for an 
ion (string) completely cooled to the ground state of the trap. For Raman 
transitions, rji can be made quite small (~ 10^^) by using co-propagating 
beams such that the Rabi frequency is given just by Q (the same reasoning 
holds for microwaves). Thus in practice, only single-qubit operations on optical 
qubits suffer from finite temperature effects. 

The speed of single-qubit rotations is often limited by the acousto-optical 
modulator used to control the light field: during switching the light field, 
phase and amplitude chirps appear which can spoil gate fidelities considerably. 
Further speed limitations are set by the amount of available laser power and 
more fundamentally by transitions close by, e.g. due to other Zeeman levels. 
Excitations on unwanted transitions due to strong laser fields transfer the 
population out of the computational subspace and thus spoil the fidelity. 



2. 5. 1 Individual addressing of ion qubits 

: Individual addressing of qubits is particularly difficult if multi and single- 
qubit gates have to be carried out in the same trap: reasonable gate times 
are only possible for high trap frequencies resulting in small ion-ion spacings, 
which in turn require very well-shaped laser beams to address the ions. 

The most straightforward way to achieve single ion addressing is to focus the 
qubit manipulation laser strongly. In the Innsbruck experiments a waist of 
2 yum FWHM reliably distinguishes between ions spaced by 5 fim. Including 
deviations from a Gaussian beam shape, a 1000-fold reduction of the light 
intensity at the position of the adjacent ions as compared to the addressed 
ion can be achieved. This corresponds to an unwanted Rabi frequency of 0.03 
on the adjacent ions with respect to the addressed ion (addressing error e = 
^cighbormg _ g_Q3^^ \ complication arises from the fact that one would like to 

^ ^addressed 

irrdiate the ion string with the laser beam perpendicular to the trap axis to 
resolve the ion positions. This in turn leads to a strongly reduced coupling of 
the laser field to the ion motion along the axial direction. In Innsbruck, an 
angle of about 68° between the laser beam and the trap symmetry axis was 
choosen as a compromise between the two diametrial requirements. 

To fulfill the stringent requirements for quantum computation, one could use 
only every second ion thus increasing the qubit-qubit distance. However, with 
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Fig. 8. Evolution of the Bloch vector on the Bloch sphere for a bit flip with 
the pulse sequence {tt /2,'k)Rz{'k)BP {it (left-hand side, see text for an 
explanation). The right-hand side displays the state evolution for a Rabi fre- 
quency corresponding to 0.3 of the one on the left-hand side (pulse sequence 
i?<^(0.37r/2, 7r)ii^(0.327r)ii'^(0.37r/2, 0)). 

an increasing number of ions, the Lamb-Dicke factor tends to become smaller 
which in turn results in slower two-qubit operations. Furthermore, the more 
complicated normal mode spectrum might reduce the obtainable fidelity. 



A more favorable method to reduce the effect of addressing imperfections is to 
use composite pulses discussed also in Sec. 12.7.31 here a single pulse is split into 
several smaller pulses. As an example, we consider here the realization of an 
i?'"(7r, 0) operation with the pulse sequence {ti /2,ti)Rz{'k)R^ {t[ /2,Q) (see 
Fig. E]). The first pulse (i?(7r/2,0)) moves the Bloch vector of the addressed 
qubit into the equatorial plane, while the Bloch vectors of the unaddressed 
qubit i is rotated by the addressing error angle ej7r/2. The AC-Stark pulse 
Rz{t^) (e.g. implemented with the laser tuned off resonance) rotates the phase 
of the addressed qubit by vr, while the phases of all other qubits are rotated 
only by e^TT. The last i?(7r/2,7r) pulse finalizes then the rotation on the ad- 
dressed qubit, while on the other qubits the effect of the first pulse is undone 
to a large extent for sufficiently small ef. Fig. [8]d illustrates this behavior on 
the Bloch sphere for a rather large addressing error of 0.3. We note that if the 
addressing imperfection of t he operations i s well- known, much better results 
can be achieved. For instance iHaljan et al.l ( j2005bl ) choose the length of a spa- 
tially inhomogeneous Rz pulse such that two qubits acquire a phase difference 
of vr (ILed . 120061 ). In this way one can afford addressing error s e close to unity. 
Similar techniques were also used by iKielpinski et al.l (120011 ). 



This composite pulse technique uses the fact that single-qubit operations com- 
mute on different ions. This is not the case for sideband transitions as multiple 
ions interact simultaneously with the a vibrational mode. However, it can be 
shown that by splitting sideband pulses into a sufficiently large number of 
pulses (~ 5) errors i n the popula tion of the vibrational mode are undone 
before they get large (iHansell . l2003l ). Thus, addressing errors can also be sup- 
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pressed in this case wit h composite pulses. A r elated technique is address- 



ing in frequency spa ce (jStaanum and Drewseru . l2002l : ISchrader et all |2004J : 
Haljan et al.l . l2005bl ): here an inhomogeneous laser or a magnetic field gra- 
dient induces different transition frequencies for each qubit. Thus the laser 
predominantly interacts with the qubit whose transition frequency matches 
the laser frequency. The disadvantage, however, is that in this case one must 
keep track of the phases of all qubits individually. In addition, it is necessary 
to control the strength of the gradient field sufficiently well to avoid dephas- 
ing of the qubits. This can be a serious problem when the qubit register is 
large and some of the qubits experience a very fast phase evolution due to the 
gradient field. 



There are various oth er possibilities to achieve effective addressing (see e.g. 
Wineland et al.l (119981 )): One particularly useful trick consists in changing the 
distance of the ions between the quant u m operations by alte ring the trap 
stiffness (see e.g. Refs. lRowe et al.l f l200lh : iReichle et al.l (l2006bh l In this way 
the relative phase of the operations on the individual ions can be changed 
allowing one to address even individual groups of ions for non-local operations. 
In addition, it was demonstrated that a two-ion string can be placed in such 
a way in the trap that the ions have a different micromotion amplitudes. 
This in turn leads to di fferent transition str e ngths and t hus to single particle 



addressing capabilities (jMonroe et al.l . Il999l : iKind . Il999l ) 



2.6 Two- qubit gates 



One route to achieve a universal set of gates (lDeutschl . ll989l ). is to complement 
single-qubit operations with two-qubit gate operations. These operations are 
one of the most important ingredients of a quantum computer as they provide 
the possibility to entangle two qubits. In combination with sin gle-qubit opera- 



tions they allow for implementation of any unitary operation (iBarenco et al. 
19951 ). 



In many implementations of quantum algorithms with trapped ions, the fi- 
delity of the whole algorithm was limited by the fidelity with which the 
two-qubit operations were implemented. Thus, currently the implementation 
of high-fidelity entangling gates is of crucial importance. In the following, 
we will discuss the Cirac-Zoller gate (Sec. 12.6.21) . the M0lmer-S0rensen gate 
(Sec. I2.6.3P and the so-called geometric phase gate (Sec. 12.6.41) . before we 
briefiy summarize various additional proposals for two-qubit gates (Sec. 12.6.3]) . 
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2. 6. 1 Motion of ion strings 



The interaction between the ionic qubits can be mediated by motional degrees 
of freedom that serve as a quantum bus for distributing quantum information 
between the ions. Therefore, we will discuss firs t the manipulat i on of t he mo- 
tion of single ions and ion strings (see also Ref. iLeibfried et al.l (l2003al )). The 
basic level scheme representing a single ion coupled to a single motional mode 
is depicted in Fig. [H A laser field resonant with the carrier transition of fre- 
quency Wqubit drives transitions with Rabi frequency VLeS where the motional 
state is not changed (see Eq. [TT]). If the laser field, however, is detuned by the 
trap frequency towards higher energies, the so-called blue sideband is excited 
(see Eq. [5]) and the operation 



R+{9,^) = exp 



i^ (e^'^a+a^ + e"'^ 



cr a 



(12) 



is carried out. Here cr^ are the atomic flip operators which act on the electronic 
quantum state of an ion by inducing transitions from the \g) to |e) state and 
vice versa (notation: = \e){g\). The operators a and denote the annihi- 
lation and creation of a phonon at the trap frequency u, i.e. they act on the 
motional quantum state. As in Eq. [10], the parameter 6 depends on the strength 
and the duration of the applied pulse, and (p denotes the relative phase be- 
tween the optical field and the atomic polarization. Importantly, the electronic 
and motional degree of freedom change simultaneously. Similarly for the op- 
posite detuning the red sideband can be excited. Fig. [9] shows a spectrum of a 
single trapped ^°Ca"^ ion near the qubit transition. At a detuning of ± 1 MHz, 
the red and the blue axial sidebands appear, respectively. In addition, radial 
sidebands (detuning ~ ±5 MHz) and higher order sidebands are visible. For 
the applied laser power and excitation time, the carrier transition is strongly 
saturated while the sidebands are only weakly saturated. This indicates that 
the sideband transitions are weaker than the carrier transition as it is expected 
from Eq. Oand Eq. [71 For the experiment shown in Fig. [9l the single ^°Ca+-ion 
was cooled to a temperature of {nx,ny,nz) ~ (3,3,16) via Doppler cooling, 
while the Lamb-Dicke factors were {r}x,rjy,rjz) ~ (0.01,0.01,0.08) (in this par- 
ticular experiment the angle of the laser beam was given by (49.2, 49.2, 22.5)). 



For quantum logic experiments, however, the case of multiple ions is more 
interesting. Since the Coulomb interaction couples the i ons strongly together , 
it is useful to find the normal modes of the ion string (ISteand . Il997l : IJamesl . 



19981 : ISasura and Buzekl . [2002). Taking a three ion string as an example and 



concentrating on the axial direction, we find three normal modes: center-of- 
mass, breathing (or stretch) and an additional axial mode. Fig. [TD] illustrates 
the motion of each ion associated with each mode. Denoting the center-of-mass 
modes frequency ui, the breathing mode has a frequency of t02 = uji and 
the third axial mode has a frequency of cus = i/29/5 cui. Fig. [TOl indicates the 
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Fig. 9. Spectrum of a single trapped ^^Ca^ ion cooled close to the Doppler limit 
(see text). On resonance (A = 0) the strong carrier transition appears as well 
as for positive (negative) detunings A blue (red) sidebands are visible. The mo- 
tional frequencies of the single ion can be deduced from the spectrum to amount to 
(j^ ~ 27r X 1 MHz (axial frequency) and ujr^ ~ 27r x 5 MHz (radial frequencies). In 
addition, higher order modes are visible at A = mtOz i ratOr-))- 



Center-of-mass: 
Breathing: 
Third axial: 



CO 



Fig. 10. Normal modes of a three-ion crystal along the axial direction with motional 
frequencies a;,. 

relative motions of the ions. The strength with which each mode couples to the 
motion is characterized by the eigenvectors (1, 1,1)/ \/3 for the center-of-mass, 
(1, 0, — 1)/\/2 for the breathing and (1,— 2,1)/a/6 for the third axial mode 
(jjamesl . ll998l ). This means that the center ion does not couple to the breathing 
mode and that the left and right ions couple with opposite phase factors to it. 
Fig. [TT] illustrates that the breathing mode cannot be efficiently excited when 
the center ion is addressed (see Sec. 12.5.11 for addressing of individual ions). 
For further details of the n orma l modes of ion strings (including larger ion 
strings), we refer to IjamesI fll998h . 



While measuring the spectra of ion strings needs only reasonably cold temper- 
atures (e.g. ions cooled to the Lamb-Dicke limit), coherent operations on the 
sideband need usually ground state cooling of this particular mode. The reason 
can be seen already in Eq. [5] which predicts that the sideband Rabi frequency 
strongly depends on the vibrational quantum number of the corresponding 
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Laser detuning at 729 nm (MHz) 



Fig. 11. Excitation spectra of a three- ion string. In Fig. [TTh). the exciting laser is 
addressed to the center ion, while in Fig. Illb) the laser is addressed to the left ion. 
Visible are excitations of the motional modes as well as of the carrier transition. 
If the center ion is addressed Fig. [TTk). the breathing mode can not be excited 
efficiently. The residual excitation of the breathing mode is most likely due to laser 
light interacting with the outer ions (see Sec. I2.5.T] on imperfect addressing). 

motional mode and thus excitation on the sideband is incoherent for a finite 
temperature of this motional mode. On the other hand, this strong sensitivity 
of the sideband excitations to the motion will be the key to couple the ion 
qubits. To allow for coherent sideband operations, the motion of the mode in 
question is cooled to the ground state, although one could also imagine using 
any well-defined motional quantum number. Fin ally, we extend Eg. [5] to take 
all motional modes into account (c.f. Eq. [TT] and IWineland et al.l (Il998) ): 



^ts = ^(V^bus^bus) 11(1 



(13) 



Here is the effective Rabi frequency on the blue sideband, fl the ideal 
carrier Rabi frequency and nbus denotes the quantum number of the bus mode 
(i.e. the mode which will be used to couple the qubits), while r/bus labels 
its corresponding Lamb-Dicke factor. The index m runs here over all modes 
except the bus mode. Since the influence of these modes on the sideband Rabi 
frequency Q'^q is strongly reduced, they are often termed spectator modes. 



As di scussed in Refs. ( iNeuhauser et al.l . ll978l : IWineland and Itand . ll979l : lMarzoli et al. 



19941 ). ground state cooling in ion traps can be achieved by driving the red side- 
band (c.f. Fig. [1]). Each absorption of a photon on the red sideband followed 
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Fig. 12. Rabi oscillation on the blue sideband of the center-of-mass mode. The data 
were taken on a string of two ^''Ca"'" ions whose center-of-mass mode was cooled 
to the ground state. Only one of the ions was addressed. The population oscillates 
between the \S, 0) and the \D, 1) state of the addressed ion. 

by a spontaneous en iission on the carrier transition takes out one motional 
quantum of energy (IStenholml . Il986l : lEschner et al.l . l2003l ). It usually takes 
too long to wait for the spontaneous decay of the upper level. Therefore in 
most experiments, the lifetime of the upper state is shortened by quenching 
the state with a laser connecting it to a fast decaying state. 



Experimentally, sideband cooling to the motional ground st ate was demon- 



strat ed first by the NIST-group with a single mercu ry ion (jDiedrich et al. 
1989) and t hen l a ter with various other ion species ( Monroe et al. ( 1995bh 



Roos et al.l fll999l ): IPeik et al.l fll999h : for a review see lLeibfried et al.l fl2nn8al )). 
For ion strings sideband cooling works simil arly as for sing l e ions. The multiple 



normal modes can be cooled sequentially (iMonroe et al.l . Il995bl : iRoos et al.l . 
19991 ). However, for ion crystals, the Lamb-Dicke factors rji tend to be smaller 



due to the increased mass as compared to single ions (for a definition of t] see 
Sec. 12. 2p and thus the cooling process is slowed down. 

Once the motion of a particular mode of the ion string is cooled to the ground 
state, irradiation on the respective blue sideband leads to Rabi oscillations 
(assuming the electronic degree of freedom of the ion is also in the ground 
state). Fig. [12] shows such a Rabi oscillation on the blue sideband of the 
center-of-mass mode of a string of two ^°Ca"*" ions. 



Another promising route to cool ions below th e Doppler limit is based on elec 



tromagnetically induced transparency (KIT) (IRoos et al.l . l2000l : iMorigi et al. 



20001 ) . Here two interfering laser fields are used to shape the atomic absorption 
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Fig. 13. Graphical representation of the three steps to perform a phase gate between 
two ions with an electronic ground state 1^), an excited state |e) and an auxilary 
state I a). 

spectrum in such a way that sharp features appear in it which can be used to 
suppress heating effects on the carrier and blue sideband transitions while still 
maintaining absorption on the red sideband transitions. In general, EIT is a 
very versatile tool as the free parameters of the two laser intensities and de- 
tunings from the atomic resonance allow one to adjust the atomic absorption 
spectrum to the actual needs (equilibrium temperature and cooling speed). 
Furthermore, iRoos et al. achieved simultaneous cooling of up to three 

modes close to the ground state with a single setting. 



2.6.2 Cirac-Zoller gate 



Conceptually, the Cirac-Zoller phase gate is the si mplest of the two-qubi t 
gates presented here and is therefore discussed first (ICirac and ZoUeii . Il995l ) . 
It requires single ion addressing and ground-state cooh ng of the bus mode. On 
the other hand, for t he M0lmer-S0rensen t ype ga tes (ISorensen and Molmerl . 



19991 : iMilburn et all . 120001 : iLeibfried et al.l . l2003bl ) both ions are illuminated 



simultaneously and the ion string have to be cooled only into the Lamb-Dicke 
limit. The geometric phase gate (ILeibfried et al.l . l2003bl ) is closely related to 
the M0lmer-S0rensen gate. However, it will be presented separately as its 
implementation looks different. 

Cirac and Zoller proposed the follo wing procedure to perfo rm a two-qubit gate 
between two ions in an ion string (jCirac and Zolleii . 119951 ): First the quantum 
information of the first qubit is transferred to the motional degree of freedom of 
a mode (the bus mode) of the ion crystal. Importantly, the resulting motional 
state affects not only the addressed ion itself but the ion string as a whole. 
Therefore on a second ion, operations conditioned on the motional state of the 
ion string can be carried out. Finally the quantum information of the motion 
is mapped back onto the first ion of the string. 



The individual steps of the Cirac-Zoller phase gate are as follows: A laser pulse 
directed to the first ion with length 9 = tt and frequency cjqaut — i^trap (red 
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sideband) moves all population present in the |0, e)-state to |1, g)-state (see left 
panel in Fig. [T3l) . However, if the first ion is in the |0, (y')-state, no state transfer 
happens. Thus, we have effectively mapped the quantum information from the 
electronic degree of freedom to the motion. Note that the coupling strength on 
the sideband depends strongly on the vibrational excitation. Therefore, this 
procedure works only if one knows the vibrational state such that one can 
properly adjust the laser intensity to achieve a vr pulse. With the quantum 
information of the first ion in the motion, we can exploit the common motion 
of the two ions and perform a c onditional operation on a the second ion. We 



follow here the original proposal (ICirac and Zollerl . ll995l ) and use a 27r rotation 



between |1, g) and a third auxiliary state |0, a) on the red sideband (see center 
panel in Fig. [13]). Importantly, only the \l,g) state is coupled to another level; 
for the states |0, e), |1, e) and |0, g), however, the coupling vanishes as there are 
no levels present with the appropriate energy. Effectively, we have therefore 
performed the following operation: 
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|0,e)^ |0,e) 

|l,e)^ |l,e) (14) 

10,^?)^ \0,g) 
\l,g) -\l,g) 

This implements a phase shift of the second ion conditioned on the motional 
state. Finally, another tt pulse on the red sideband addressed to the first ion 
maps the quantum information present in the motion back onto the first ion 
and in total a phase gate between the two ions is performed. To turn this 
phase gate into a controlled-NOT operation, one can apply an {ji /2,ipi) 
pulse {i = 1,2) to the target ion on the carrier transition (see Eq. [TOj) before 
and after the phase gate. Choosing a particular phase relation (^2 — ^1 either an 
ordinary controlled-NOT gate or a zero controlled-NOT operation is obtained. 

The heart of this procedure, i.e. the dynamics presented in Eg . [T^ was demon - 



strated by the NIST group with a single beryllium ion (iMonroe et al.l . Il995al ). 
The Innsbruck group impler nented the complete protoco l with individual ad- 
dressing of two calcium ions fISchmidt-Kaler et all . l2003cl Jbl) . Sec. O describes 
both experiments. 



Jonathan et al. fl2000h present a generalization of this gate where they propose 



to drive the system so strongly that due to AC-Stark shifts the eigenstates 
of different vibrational quantum numbers (the dressed states) get the same 
energy. Thus a motional dynamics can be achieved just with carrier pulses of 
appropriate phase. Essentially, this proposal trades faster gate speeds against 
more sensitivity to laser intensity noise as the energy of the dressed states 
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Fig. 14. Energy-level diagram of two trapped ions illustrating the principle of the 
S0rensen and M0lmer gate. The bus mode is populated with n phonons. Two laser 
beams tuned close to the blue and red sideband, respectively, drive the system via 
the dashed virtual levels between the \n,gg) and \n,ee) state. A similar process 
takes place if the ion string is either in the |n, eg) or in the |n, ge) state. 

depends strongly on the laser power. 



2.6.3 M0lmer-S0rensen gate 



Another possibility to implement a t wo-qubit gate is to use laser radiation 



tuned close to the motional s i debands fISeirensen and Meilmerl . Il999l ; iMeilmer and Sgirensen 
1999 : S0rensen and Maimer . 2000l : Solano et al. . 19991 ). The basic principle is 
as follows: both ions are irradiated with a bichromatic laser field with fre- 
quencies Uq ± (wqubit + S), tuned close to the red and the blue sideband of 
a collective mode, respectively (see Fig. [H]). The two frequencies sum up to 
twice the qubit frequency cUqubit, each laser field itself, however, is not reso- 
nant to any level. Thus both ions can change their state only collectively and 
choosing an interaction time appropriately, the dynamics 



|ee)^(|ee) +i|gg))/V2 
|eg)^(|eg) + i|ge))/V2 
|ge)^(|ge) + i|eg))/V2 

|gg)^(|gg) + i|ee))/y2 (15) 

is achieved. To see that this gate leads to a universal set of gates, we introduce 
the new basis |±)j = (|e)i ± \g)i)/\/2. |±)|±) are eigenstates of the unitary 
operation described by Eq. [TSland transform as |++) — > 1++), | + ~) ^ "^IH — ); 

I — h) \ — h), I ) I ) by the action of the gate, where we omitted 

a global phase factor of e"*'^/'^ on the right hand sides. This transformation is 
a conditional phase gate up to single-qubit phase-shifts and is known to be 
universal together with single-qubit operations. 
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The M0lmer-S0rensen gate has the particular feature that it does not require 
individual addressing of the ions and that it does not fail completely if the ion 
string was not cooled to the ground state. As will be described i n Sec. 14.21 the 



NIST group entangl ed up to four ^Be"*" ion s with this gate type (ISackett et al. 



2OOOI ). Furthermore, iHaljan et al.l (l2005bl ) created all four Bell states and im- 
plemented Gro yer's seach algor it hm with ^^^Cd"*" i ons us ir ig this gate opera tion 



see also Refs. iLee et all fl2005f ): iBrickman et al.l (120051 ): iBrickmanI (120071 )) 



Very recently, the application of a M0lmer-S0rensen gate t o optical qu bits 



excited on dipole-forbidden transitions has been analyzed (iRood . l2008l ). It 
was shown that fast, high-fidelity gate operations are achievable by smoothly 
switching on and off the bichromatic laser fields inducing the gate action. The 
experimental implementation of this te chnique has resulted in the creation of 
Bell states with a fidelity of 99.3(1)% (iBenhelm et al.l . l2008bh . 



2.6.4 Geometric phase gate 

The so-called geometric phase gate uses also two laser fields irradiating mul- 
tiple ions at the same time. An interesting feature of this gate is that ideally 
during the gate operation the electronic state of the ions is not touched. Only 
a force dependent on the electronic states is ap plied such that f or the various 



electronic states different phases are acq uired (iMilburn et al.l. 120001) . In the 



scheme implemented by the NIST group (iLeibfried et al.l . l2003bl ) two non-co- 
propagating laser fields create a standing wave (see Fig. [T5|) . The difference of 
the two laser frequencies is tuned closely to one of the axial frequencies which 
leads to a walking wave. Thus each ion experiences a periodic AC-Stark shift 
and a force depending on the slope of the spatial variation of the AC-Stark 
shift. Most importantly the size and even the direction of the force can depend 
on the electronic state of the ion. Choosing the distance between the ions such 
that each ion experiences the same phase of the walking wave for a given time 
(see Fig. [T5|) . ions are pushed in the same direction if they have the same 
electronic state. In this way the breathing mode cannot be excited. However, 
if the ions are in different electronic states the forces on the ions are not the 
same and the breathing mode can be excited. The detuning of the drive from 
the motional mode is chosen such that the phase between the motion and the 
drive changes its sign after half the gate time. In this way, the ion string is 
driven back to the original motional state after the full gate time. This en- 
sures that the motion is disentangled from the electronic state after the gate 
operation. The intermediate energy increase as compared to situations where 
the motion is not excited leads to the desired phase factor and the ions have 
picked up a phase factor that depends non-linearly on the internal states of 
both ions. 
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Fig. 15. Force on two ions in a standing wave of two laser fields, a): both ions spin 
up. Both ions experience the same force b) one ion spin down, the other spin up. 
Both ions experience opposing forces. Tuning the frequency difference of the two 
laser fields close to the breathing mode frequency, in a) the motion of the ion string 
cannot be excited efficiently, while in b) the breathing mode is excited. 



Leibfried et al.l (l2003bl ) implemented this gate with a fidelity of 0.97, limited 
mainly by the spontaneous decay from the P manifold of the Be"*" ions during 
the gate operation. This exceptional fidelity was reached because the gate 
avoids a number of imperfections in the first place. We mention here the 
absence of off-resonant excitations on the carrier transition. Furthermore, the 
gate execution time of about 10 /is is two to three orders of magnitude faster 
than the corresponding coherence time of the qubits and thus decoherence 
effects are small. To obtain a differential force on the ions, the NIST group used 
two laser beams detuned blue from the 6*1/2 P1/2 transition by 2tt x 82 GHz 
(fine structure splitting of the P state is 200 GHz). For the ratio of the forces 
on the 1 1) = |F = 2,mj = —2) and the 1 1) = |F = 1,mf = —1) state 
this yields: Fi=-2Fp Additional Stark shifts can be efficiently suppressed by 
choosing almost perpe ndicular a nd linear polari z ation s for the laser beams 



(jWineland et al.l . l2003l ). Finally, iLeibfried et al.l (120071 ) discuss a version of 



this gate where the laser intensities impinging on the ions are controlled by 
transporting the ion crystal through the laser fields. Controlling the interaction 
of the ions by transport offers also the possibility to use spat ially modulated 



magn etic fields to create the state dependent oscillating force (ILeibfried et al. 



20071 ). 



These three gate types have different strengths and weaknesses. Both, the 
geometric phase gate and the M0lmer-S0rensen gate do not need single ion 
addressing. While this is often advantageous, it has the inconvenience that 
it is not straightforward to carry out these gates between specific ions in a 
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string. Using segmented traps for moving and splitting ion strings (see Sec. [6]) 
resolves this issue. Another route to introduce ion specific gates is to hide 
certain ions with single ion operations (see Sec. 14. 5p such that these ions are 
not affected by the multi-qubit gate operation. The Cirac-Zoller gate on the 
other hand demands single ion addressing but allows for a straightforward 
implementation of quantum algorithms. 

Another aspect is speed. The geometric phase gate can be executed quite fast 
as the laser can be tuned such that off-resonant transitions are quite unlikely. 
The Cirac-Zoller and M0lmer-S0rensen gates on the other hand require a laser 
tuned to or tuned close to a sideband transition. This automatically implies 
that the laser is detuned only by (approximately) the trap frequency from the 
strong carrier transition. Thus it seems that the gate speed has to be much 
slower than a trap frequency (see Sec. 12.21) . However, for special temporal and 
spectral arrangements of the laser field, it is possible to suppress the spectral 
contribution on the carrier transition so that gate times close to the trap 
period seem feasible. 



Finally, it is important whether the gate works efficiently with magnetic field 
insensitive transitions. Qubits encoded in superpositions of levels connected 
via a magnetic field insensitive transition provide very long coherence times 
of many seconds (Sec. 13.1.21) . In its originally proposed form, the geometric 
phase gate is very inefficient on magnetic field insensitive transitions because 
hyperfine states w i th a s imilar magnetic moment appear to experience similar 
AC-shifts (iLangeij - liooeh . However, this could be circumvented by using a pair 
of laser beams tuned close (as compared to the qubit difference fr equency) to a 



spect rally narrow transition to induce the state-dependent force (lAolita et al. 



20071 ). In this way one induces a spectral rather then a polarization dependent 



fore. Another option is to recode the qubits for the conditional phase gate 
from a magnetic insensitive coding to a different coding (ILangerl . |2006| ). The 
M0lmer-S0rensen gate, on the ot her hand, has been a lready applied on mag- 
netic field insensitive transitions (IHaljan et al.l . l2005bl ). The Cirac-Zoller gate 
(either with composite pulses described in Sec. 12.7.31 or a third magnetic field 
insensitive transition as available for instance with the D manifold in ^^Ca"*") 
can be also directly applied in these situations. 



2.6.5 Other gate types 



Another po ssibility to employ A C-Stark shifts for a two-qubit gate was demon- 
strated by iBrune et al.l (Il994l ) with Rydberg atoms passing through a mi- 
crowave cavity. We s umm arize here the ion trap version implemented by 



Schmidt-Kaler et al.l (120041 ). In these experiments, a laser was tuned close 



to the axial motional sideband of a single ion. The Rabi frequency of the blue 
motional sideband f2+ depends on the phonon occupation number such 
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Fig. 16. Simplified level sche me of two trapped '^^Ca^ ion s including the motional 



state of the bus mode (from ISchmidt-Kaler et al.l (|2004l )). The virtual levels are 



represented by dashed lines. A monochromatic laser beam can only induce transi- 
tions between the \ge) and the \eg) state (solid arrows). In addition, the \gg) state 
acquires an additional phase factor due to an AC Stark effect (dashed arrows). 

that the resulting AC-Stark shift of the electronic levels due to the sideband 
resonance depends on the motional state: 

^E = h% = h'^{n-, + l) , (16) 

where we used Eq. [51 Rotating the frame by exp(^^t) removes the phase 
evolution of the two electronic levels in motional ground state = 0. The 
relative phase of the two electronic states in the first excited state = 1, 



however, evolves as = Thus, choosing the interaction time t = 

the phase gate operation diag(l, 1, is implemented, where we used the 

basis {|e,0), \g,0), |e, 1), \g, 1)}. 



This gate can also be generalized to multiple ions (ISchmidt-Kaler et al.l . 120041 ) . 
Here all ions simultaneously interact with the laser beam. The situation is 
very reminiscent of the M0lmer-S0rensen gate: taking as an example a two- 
ion crystal, only transitions between the eigenstates with the same number of 
excited ions can be induced, e.g. between the \ge) and the \eg) (see Fig. [T6|) . 
All other basis states acquire only a phase factor which can be corrected for 
with single-qubit operations. Using this method starting from the \eg) state, 
the Bell state {\ge) + \eg))/\/2 was generated with a fidelity close to 0.9 by 
applying a pulse corresponding to a n/2 dynamics with trapped ^"Ca"'" ions 
(H. Haffner et ai, Innsbruck, unpublished). 

Another gate type uses the dependence of the carrier Rabi freq uency from the 
motional state (see Eq. [11]). Proposed bv lMonroe et al.l (19971). it was imple- 
mented later by the same group on a single ^Be"*" ion (iDeMarco et al.l . |2002| ). 
Essentially, a laser pulse on the carrier transition is applied, with the Lamb- 
Dicke parameter chosen such that for one motional state an even number of 
Rabi oscillations occurs, while for the other motional state an odd number of 
Rabi cycles occurs. This way the phase of an electronic state is only flipped in 
the latter case, while it remains unaffected in the former case. This implemen- 
tation of a controIIed-NOT operation requires, however, a reasonably large 
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and tunable Lamb-Dicke factor of the bus mode r]b (the gate time is quadratic 
in l/r]b), while either all other motional modes are cooled to the ground state 
or their Lamb-Dicke factors are much smaller than rjh- 

Finally, there have been quite a number of other promising gate proposals 
which are not implemented as of yet. Some of them rely — as some of the 
already prese nted gates — on state dependent AC-Stark shifts, albeit i n the 
static regime dCirac and Zolleil . I2OOOI : ISteand . l2004l : IStaanum et aP . l2004h . The 
general idea here is that an inhomogeneous laser beam changes the distance 
between the ions depending on their electronic state. Simulations show that 
the fidelities can be quite high and that the experimental requirements are not 
very demanding. For instance, only moderately low temperatures are required. 
Alternatively, the sta te dependent potential can be c reated with strong mag- 
netic field gradients (IMintert and Wunderlichl . 1200 ll ). Creation of sufficently 
strong field gradients might be eased by implementing these gates in micro- 
fabricated ion traps where smaller ion-surface distances and microstructured 
current carriers would facilitate the generat i on of the required static or dy- 
namic field gradients (ILeibfried et al.l . 120071 : IChiaverini and W. E. Lybargerl . 
20081 : lOspelkaus et alLboOSh . 



Another gate class, which recently received muc h attention, uses short bu t 



strong laser pulses to kick the ion string strongly (iGarcfa-RipoU et al.l . |2003| ) . 
Most interestingly, the corresponding gate operation times can be shorter 
than even one trap period. Choosing proper phases and amplitudes of the 
pulses, the ions are kicked in such a way that they acquire a state dependent 
phase due to their motion. Anot her gate variant use s a co ntinuous irradia- 
tion w'ith fast phase modulatio ns (IGarcia-Ripoll et al.l . |2005| ). Based on these 
ideas, iDuan and Kimbld (120041 ) propose quantum computation in an array of 
trapped ions where the fast gates are used to induce a next neighbor interac- 
tion. 



2.7 Apparative requirements 



Experimental quantum computation requires exceedingly long coherence times 
as well as an exquisite control over the qubits. The former is achieved by decou- 
pling the qubits from the environment, while the latter requires a well-defined 
and switchable interaction with the environment. In ion trap quantum com- 
puting, the qubits are held in free space with electromagnetic forces which 
hardly couple to the ion's internal degree of freedom. The control is achieved 
with strong laser fields such that during their interaction they still can be 
treated classically and entanglement between the laser field and ions is negli- 
gible. Thus the seemingly contradictory requirements — decoupling from the 
environment and controlled interactions — can be fulfilled. 
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Fig. 17. Front (le f t) and side (right) view of the Innsbruck trap (from 
Schmidt-Kaler et al. ( 2003bl )). It consists of four blades (light gray) and two end 



caps (dark gray). A radio- frequency drive is applied to two opposing blades while 
the other blades are held at ground. This provides confinement perpendicular to the 
trap axis. The two end caps held at a positive DC-potential U^nd prevent the ions 
to escape along the axis. 

For ion-trap quantum computing usually so-called linear Paul traps are used 



wher e the ions can form a linear ion chain (iRaizen et al.l . Il992bl : iDrees and Paul! . 



19641 ). In these devices, a radio- frequency potential is applied to two electrodes 
which are parallel to the axis of the trap (see Fig. [T7j). These electrodes create 
an oscillating two-dimensional quadrupole potential that is translation invari- 
ant along the trap axis. If the frequency of the RF field is sufficiently large, 
the ions experience an effective restoring force to the center axis. Additionally, 
static electric fields confine the ions along the trap axis. The trap resides in 
an ultra-high vacuum vessel to reduce collisions with residual molecules and 
atoms as much as possible. For typical experiments, it is sufficient to describe 
the time dependent confining forces of the ion trap as if they resulted from a 
static three-dimensional harmonic potential. If the confinement perpendicular 
to the trap axis (radial direction) is much larger than the confinement along 
the trap axis, cold ions form a linear chain. Typical trap frequencies for the 
radial frequencies are between 4 and 10 MHz, while axial frequencies range 
mostly between 0.5 and 5 MHz. 

For the control of electronic and motional states of the ions, lasers with high 
frequency and intensity stability are used. Acousto-optical modulators allow 
one to control both the frequency and the intensity. Typically, the lasers are 
either referenced to an ultra-stable cavity or to a molecular transition. Ions 
tend to have higher energy splittings as compared to atoms. Therefore often 
lasers emitting in the ultraviolet are required. In the next paragraphs, we 
briefly describe the experimental setups used in the NIST and the Innsbruck 
experiments. 



2.7.1 The NIST setup 



Details of the N IST setup can be found in the PhD theses of iLangen (|2006[ ). 



Kielpinskil ( 120011 ) and lKing] (119991 ). The NIST group uses ^Be"*" ions whose level 



scheme is depicted in Fig. [181 Qubits are encoded in the hyperfine manifold of 
the Si/2 electronic ground state split by 1.25 GHz. For most of the experiments 



34 



P state 



m, = 1/2 



= -3/2 

200 GHz 



240 MHz 



F = 1 



F = 2 



313 nn 



S state 



F = 1 



1.25 GHz 



-1 m - 2 

= 2 F p ^ 2 



Fig. 18. Energy level scheme of a ^Be"*" ion (nuclear spin: 3/2). The hyperfine split- 
ting of the P3/2 state is smaller than 1 MHz and not shown. 



2, rriF 



-2) 



discussed here, the qubit was encoded in the \F 
lyTUp = — l)-transition as the levels are easily prepared and distinguished 
from each other. Doppler cooling at 313 nm (laser power ~ /iW at a waist 
of ~ 25 /im) and optical pumping with a~ polarized light on the \Si/2,F = 



2, nip 



-2) 



IP' 



3/2, 



3, mp 



-3) transition along the quantization 
axis given by a weak magnetic field (typically i? ~ 1 mT) initializes the ions 
in the \F = 2,mp = —2) state. Finally, pulsed resolved sideband cooling 
on the \F = 2,mF = —2) < — > \F = 2,mp = —1) transition is used to 
prepare the ions' motio n with a probability of about 0.99 in the ground state 
(jWineland et al.l . 119981 ). before the quantum information is manipulated. 



The qubit states are coupled with a Raman transition via the P manifold. 
The NIST group generates the necessary Raman beams from the same laser 
source such that the relative phase of the beams is well-defined. Therefore the 
dominant decoherence source is dephasing due to magnetic field fluctuations 
(e.g. due to the mains supply at 60 Hz) causing a qubit lifetime on the order 
of a few milliseconds. However, recently the NIST group used a magnetic field 
insensitive transition and measured coherence times on the order of a few 
seconds (see Sec. 13.1.21) . 

Read-out is performed again on the cycling transition \Si/2,F = 2,mp = 
—2) i — > \P:i/2,F = 3,mp = —3) with a~ polarized light. This light does not 
couple efficiently to the \Si/2,F = l,mp = —1) state as it is detuned from 
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any possible transition by a little bit more than 1 GHz. In order to avoid 
pumping into the bright state by the still present off-resonant excitations, 
the |S'i/2,-F = l,mp = —1) population is transferred with two vr pulses to 
the 15*1/2, -F = iymp = +1) level. In this state, absorption of a single off- 
resonant photon cannot lead to a population of the l>S'i /2, F = 2 , mp = —2). 
Thus the ion remains dark. The theoretical analysis by iLangerl (120061 ) shows 
that with this method detection errors can be kept smaller than 10~^. At the 
moment, due to stray light background and imperfections in the vr transfer 
pulses, typically dete ction efficiencies of about 99% are attained in the NIST 
setups (ILangerl . |2006| ) . 



The small atomic mass of beryllium allows for high trap frequencies and com- 
paratively large Lamb-Dicke factors. Both factors alleviate a strong coupling 
to the motional degree of freedom and thus facilitate fast multi-qubit opera- 
tions. In order to initialize, manipulate and detect the quantum states, light 
sources in the ultraviolet at 313 nm are required. 

The NIST group uses various microstructured traps. Many of these traps have 
multiple trapping zones, whi ch ease the scaling of ion trap quantum com- 
puters to larger ion numbers ( Kielpinski et al.l . 2002). In addition, segmented 
traps allow for single-qubit addressing without tightly focused laser beams 
(see Sec. 15. 2p and a separate loading zone avoiding patch effects (see Sec. | 



2.7.2 The Innsbruck setup 



The experimental se tup used by Innsbruck group is described in Refs. (ISchmidt-Kaler et al. 



2003bl : iGuldd . l2003l ). As qubits, superpositions of the S1/2 ground state and 
the metastable D5/2 state of ''°Ca+ are used (see Fig. Hj). The D5/2 state has a 
lifetime r ~ 1.16 s. 



A magnetic field of 300 /xT lifts the degeneracies of the Zeeman manifolds. For 
the experiments, the entire quantum register is prepared by Doppler cooling, 
followed by sideband cooling to the motional ground state. Normally, only the 
center-of-mass mode (cc'cm ~ 1-2 MHz) is cooled to the ground state. The ions' 
electronic qubit states are initialized in the Si/2(^«j = —1/2) state by optical 
pumping with a~ light. Then each ion-qubit is individually manipulated by 



a series of laser pulses on the S = S1/2 



m,- 



=-1/2) to D = D 



5/2 



m^ 



=-1/2) 



quadrupole transition near 729 nm. In order to guarantee phase coherent ma- 
nipulation, the laser frequency is stabilized to about 50 Hz on a time scale of 
1 minute by loc king it to an ul t ra-sta ble reference cavity with a similar design 
as presented by lNotcutt et al.l (120051 ). This time scale corresponds to the typ- 
ical duration of a full set of quantum computing experiments. On time scales 
of a few seconds even a 3 Hz linewidth has been observed. The slow drift of 
the reference cavity — typically about 1 Hz/s — is monitored every few min- 



36 



utes by interrogating the qubit transition and compensated with a feedback 
loop. Thus many thousands of experiments under comparable conditions are 
feasible. 

Spectroscopy on a transition more sensitive to magnetic field fluctuations (S1/2 
(mj=-l/2) —>■ D5/2 (mj=-5/2)) is used to monitor slow magnetic field drifts 
continuously. In addition, a passive magnetic shield reduces the magnetic field 
fluctuations on timescales longer than 1 s by about a factor of 30 and by more 
than 2 orders of magnitude for frequencies higher than 10 Hz. The total field 
amplitude noise is on the order of a few nT. Thus, typically coherence times 
of about 5 ms on the S1/2 (mj=-l/2) D5/2 (mj=-5/2) transition and 15 ms 
on the Si/2 (k1j=-1/2) — > D5/2 (mj=-l/2) transition are achieved. The laser 
driving these transition is tightly focused onto individual ions in the string 
with a waist size of 2 /im (inter- ion distance ~ 5 /im). 

While the NIST group uses predominantly global addressing and state read- 
out, the Innsbruck-group uses tightly focused laser beams to address individual 
ions (see Sec. I2.5.ip . Thus in principle any quantum algorithm can be imple- 
mented straightforwardly, only limited by the decoherence time. The trade-off 
of this method, however, is that the axial trap frequency can not be increased 
too much, since then the ions move closer to each other, thwarting single ion 
addressing. A consequence of a lowered trap frequency is a reduced speed of 
the entangling operations on the sidebands. 



2.7.3 Composite pulses and optimal control 



Quantum algorithms are usually implemented by a sequence of a few funda- 
mental gates. Many of those gates can be carried out with single laser pulses. 
However, using a set of pulses sometimes offers an advantage over using sin- 
gle pulses (as demonstrated already in Sec. 12.5.11) . In NMR, the composite 
pulse technique is well known and allows for the compensati on of many sys- 
tematic effects like intensi t y mismatch and detuning errors (IFreemanl . 119971 : 
Levitt an d Fre emanl . 1 19791 : iLevittl . Il986l ). The spin echo sequence discovered 
by iHahn (1950) is such a sequence with which a constant detuning between 
the excitation and the transition can be removed to a large extent. 



In addition to the many sequences discovered and used in NMR, there are a 
few sequences which are important in the context of ion traps. We present 
here two of thos e which were used in the implementation of the Deutsch- 
Jo sza algorithm (Guide et al. . 20031 ) (see Sec. 15.11) and were described first 



bv IChilds and ChuaneJ (l2000h. The firs t sequence uses four sideband pu lses 
flChilds and Chuand . I2OOOI : ICulde et all . l2003l : ISchmidt-Kaler et all . l2003bh to 
implement a phase gate in the computational subspace {\D, 0), 15", 0), \D, 1), IS", 1)}. 
The advantage over the method laid out in Sec. 12.6.21 is that no third level is 
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Fig. 19. Evolution of the Bloch vector during the composit e phase gate in the 



\ S. 0) ^ \D, 1) (left) and 15", 1) ^ \D, 2) manifold (right) ffrom lSchmidt-Kaler et al 

jSoS)). 



required to implement the gate. To achieve the desired gate with just two lev- 
els, cleverly chosen pulse lengths and phases avoid leakage into higher phonon 
number states. Using the definitions from Eq. [121 the pulse sequence (to be 
read from right to left) is: /?+(7r/2, n/2)R+{TT/V2, 0)R+{n/2, n/2)R+{TT/V2, 0). 
Having Fig. [1] in mind, we analyze the effect of this pulse sequence on the four 
physical eigenstates {\D,0),\S,0),\D,1),\S,1)}. The \D,0) state is not af- 
fected at all and therefore \D, 0) \D, 0). Fig. [19] illustrates the evolution of 
the respective Bloch vector in the 15, 0) ^ \D, 1) and 15*, 1) ^ \D, 2) manifold, 
respectively. Because the couplings and thus the effective pulse lengths differ 
by a factor of between the two manifolds, the Bloch vector follows different 
paths. Still, it reaches always its original position when the pulse sequence is 
finished. Using Eq. [12] and Eq. [12], one can show that for the three remaining 
cases, each time a phase factor of -1 is picked up and thus the diagonal matrix 
diag(l, —1, —1, —1) is implemented. 

Similarly, three blue sideband pulses can be employed to implement a SWAP 
oper ation between an e lectronic and a motional degree of freedom of trapped 



ions flGulde et all . 120081 ): R+{7r/V2, 0)i?+(27r/V2,(^swap)i?+(7r/V2, 0), where (p, 



arccos (^cot^(7r/\/2)) . This pulse sequence was used by Guide et al.l ( 2003 ) to 



swap 



implement the Deutsch-Josza algorithm (Deutsch, 19891 ). Already these two 



pulse sequence examples demonstrate that composite pulses are a quite ver- 
satile tool. 

Especially, gradient ascent pulse engineering (GRAPE) develo ped in the con- 



text o f NMR carries the idea of composite pulses to its extreme (iKhaneja et al. 



20051 ). Here a pulse sequence thought to implement a particular unitary op- 
eration is split into many pulses. Then a special algorithm is used to vary 
amplitudes and phases of the pulses to perfect the desired unitary that is 
optimal with respect to certain quality criteria (e.g. execution time) with a 
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special algorithm (jKhaneja et all l2005l ). In addition, various boundary con- 
ditions (e.g. experimental constraints like finite pulse rise times) can be in- 
cluded in form of cost functions. Thus, significant performance improvements 
in terms of speed and reduced susceptibility to experimental imperfections can 
be achieved. Most interestingly, the sensitivity to control parameters can be 
minimized with GRAPE, too. However, note that in spite of the fact that the 
execution time can be reduced considerably, composite pulses lead usually to 
a larger total pulse area. Thus decoherence effects which scale with the pulse 
area might become relevant. One such source of decoherence is spontaneous 
emission during Raman transitions (see Sec. 13.1.31) . 

First steps in apply ing such optim a l cont rol techniques to trapped ions have 
already been taken. iTimoney et al.l (120081 ) deviced pulse sequences to robustly 
perform 7r/2 and tt rotations between two hyperfine qubits of a single ^^^Yb"*". 
Creating the states |0) + e*'''|l) and |1) from |0), this work also experimentally 
demonstrates the robustness of these pulse sequences as compared to simple 
7r/2 or TT pulses, respectivvely, with respect to intensity and detuning errors. 



Furthermore, iNebendahl et al.l (120081 ) modify a GRAPE algorithm to con- 
struct a controlled-NOT operation from a global M0lmer-S0rensen interaction 
and single qubit operations. The algorithm allows for the optimization of whole 
algorithms. Taking as an example a quantum error correction scheme for bit 
flips based three qubits and two ancilla qubits, the original length of more 
than 100 pulses was reduced to 34 pulses. 



3 Decoherence in ion trap quantum computers 



This section describes the most relevant decoherence mechani sms for ion trap 



juantum computers. For further discussions, we refer to Ref. IWineland et al. 
1l998h . 



3. 1 Sources of imperfections in ion trap quantum computers 



3.1.1 Bit- flip errors 

Bit-flip errors occur when a process transfers populations between the phys- 
ical eigenstates of the qubit. Usually, the physical eigenstates are the energy 
eigenstates of the system and bit-flips are connected to radiation or absorption 
of photons. Thus, bit-flip errors are usually caused by spontaneous emission. 
The frequency differences of hyperflne and Zeeman qubits are quite small and 
therefore the time constants for spontaneous emission are usually longer than a 



39 



X -50 

LU 





Laser detuning at 729 nm (kHz) 



50 



Fig. 20. Ramsey fringes with a free precession time of T = 100 /is (Chwalla et al, 
Innsbruck, unpublished). The data points indicate the excitation probabihty to the 
metastable D level in a single ^^Ca+ ion. 

year and thus are of no relevance if the ions are not irradiated with electromag- 
netic radiation. For optical qubits usually superpositions of the ground state 
and the metastable state D-level in earth- alkali elements are u s ed. Typical life 
times of the D-levels are about 1 second ( Barton et al. . 2000l : Kreuter et al 



20051 : iLetchumanan et al.l . |2005| ). which is long as compared to the gate time 



of less than 1 ms. Therefore, most current experiments are not yet limited 
by bit-flip errors during their free evolution. Bit-flip errors on the motional 
degree of freedom are discussed in section 13.21 



3.1.2 Dephasing 

The phase evolution of a superposition often depends on a classically well- 
defined parameter such as the magnetic field. Ignoring the time evolution of 
the classical parameter leads to dephasing of the superposition. For instance, 
a superposition consisting of two levels with differing magnetic moments ex- 
periences energy shifts due to a (fluctuating) magnetic field. Thus the phase 
of the superposition depends on the particular magnetic field history and de- 
phasing occurs. Similarly, dephasing also takes place for the ion motion if the 
trap frequency is not constant, e.g. due to voltage fluctuations on the trap 
electrodes (see section . 

A typical experimental sequence to test the phase coherence consists of two 
7r/2 pulses separated by a waiting time T (a Ramsey experiment). During the 
waiting time, any difference between the atomic transition frequency and the 
laser frequency will lead to an evolution of their relative phase. The second it/ 2 
pulse will then rotate the atomic state according to the phase difference either 
towards the excited or the ground state. This effect can be seen in Fig. [20] when 
the population of the ion oscillates with the detuning of the laser frequency 
from the atomic transition (Ramsey fringes). The contrast (sslzSIH) in i^i^g 

\ '=> / \ max+mm / 

center of this curve is in the following termed Ramsey contrast. 
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Experimentally, fluctuations of the atomic resonance frequency, e.g. typically 
due to fluctuations of the magnetic fleld and of the laser frequency, lead to 
dephasing. Sometimes it is useful to distinguish between fast fluctuations (the 
relative frequency changes during the waiting time), slow fluctuations (the 
relative frequency is constant in each experiment, but not the same in the 
next experiment (s)) and an intermediate regime. It can be shown that with 
an increasing waiting time, fast fluctuations lead to an exponential decay of 
the Ramsey contr ast, while slow fluctuat ions lead to a Gaussian decay of the 
Ramsey contrast (ISengstock et al.l . Il994j ). Furthe rmore, the s low fluctuations 



can be compensated for by spin-echo techniques (Hahn, 1950l ). 



Currently, the coherence time of most qubits is limited by magnetic fleld fluc- 
tuations to a few milliseconds. Usually Fourier components at multiples of the 
mains frequency contribute most to the magnetic fleld changes. Synchronizing 
the experiment with the phase of the mains supply counters time variations 
of the magnetic fleld effectively. 



A more generic way to reduce magnetic fleld fluctuations is shielding of am- 
bient magnetic flelds with /i- metal and/or active cancellation with a feedback 
loop. Currently, the Innsbruck group uses an aluminum/ /i metal shield which 
suppresses the magnetic fleld fluctuations by more than two orders of mag- 
nitude. However, in most cases a flnite magnetic fleld is needed to lift the 
degeneracy of the electronic levels. Therefore, a very stable magnetic fleld has 
to be created inside the magnetically shielded region without saturating the 
yU-metal. The standard procedure is to use a pair of Helmholtz coils through 
which a current is passed. As it is very difficult to stabilize currents to better 
than 10~^, in the future it might be worth while going through the trouble to 
use superconducting solenoids. Here magnetic fleld stabilities of better than 
10^^^ at 6 T have been attained. These experiments were made possible by 
choosing a particular g eometry of the sup erconducting coils to shield the ex- 



ternal magnetic flelds (jPyck et al.l . 119991 ). especially of high importance for 



small flelds. It remains to be seen to what extent these exceptional fleld sta- 
bilities can be obtained at the relatively small flelds of 1 mT as required for 
quantum computation. 



A more elegant solution to reduce dephasing due to magnetic fleld fluctuations 
is to use qubit levels having the same magnetic moment. For this especially 
ions with a hyperflne structure have interesting levels. Obvious choices for the 
qubit transitions are of the form rriF = ^ mp = 0, which experience only 
a quadratic Zeeman-effect at small magnetic flelds. However, at the magnetic 
flelds required to lift the Zeeman degeneracies, a considerable linear Zeeman- 
effect is present. Therefore, it seems advantageous to work with stronger mag- 
netic flelds where transitions with a purely quadratic Zeeman-effect can be 
found. 
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While magnetic field insensiti ve transitions were extensively explored in mi- 
crowave precision experiments (jBoUinger et alJ . ll99ll : lThompsonl . [T990h . in the 
context of quantum compu ting it has been imp lemented only ver y recently by 
the N IST group on ^Be+ flLaneer et all . l2005h. by the Oxford (^Lucas et al.l 



20071 ) and the Innsbruck groups on ^ ^Ca"*" ( Benhelm et al.l . 120071 . l2008al ) and 



by the Ann- Arbor group on Cd"*" (iHaljan et al.l . l2005bl ). To achieve a rea- 



sonable spatial selectivity, the qubits are manipulated with an optical Raman 
drive (c.f. Sec. 12.31) . However, care has to be taken to maintain the phase 
coherence between both laser fields on time scales of several seconds. Using 
co-propagating laser beams simplifies this obstacle considerably as both beams 
propagate along the same path such that effects of mirror vibrations as well 
as of air currents cancel. However, a co-propagating geometry does not allow 
for easy coupling to the motional degree of freedom. On the other hand, if 
only phase gates are used as two-qubit gates, the phase stability only has to 
be maintained during each phase gate operation. 



Langer et al.l (120051 ) demonstrated coherence times r > 10 s using a magnetic- 
field-independent hyperfine transition in ^Be"*" at a magnetic field of Bq 
0.01194 T. In these experiments, Ramsey spectroscopy was carried out on 
the \F = 2,mF = 0) < — > \F = l,mF = —1) qubit transition to measure the 
phase coherence. The optimal magnetic field was determined by measuring this 
transition frequency (~ 1 GHz) as a function of the magnetic field. The min- 
imum of the resulting parabola (second order derivative B2 — 0.305 Hz//iT^, 
Bq ~ 0.01194 T) corresponds to the desired magnetic field with the least 
magnetic field dependency. The coherence time of this qubit is limited by slow 
drifts of the ambient magnetic field within the typical measurement times of 
a few hours especially for scans with long Ramsey waiting times. These ex- 
periments demonstrate a qubit memory error rate on the order oHq/t 10~^ 
where the time scale to is set by the detection time of to = 200 yus, which is 
the longest operational time of the NIST ion-trap quantum computer. 



The Oxford gr oup carried out co herence measurements on ^'^Ca+ with a mi- 
crowave drive (ILucas et al.l . 120071 ). They used the |F = 3,mi? = 0) ^ \F = 
4, mp = 0) transition in the hyperfine ground state manifold of ^^Ca'*' at small 
magnetic fields {B ^ 0.178 mT) and observed a dephasing time of 1.2 (2) s. 
Additionally, they investigated the coherence properties in a spin echo con- 
figuration and could not detect any decay of the Ramsey contrast on time 
scales of up to 1 s. From this they deduce a spin-echo dephasing time of larger 
than 45 s. Similar results were obtained by the Innsbruck group in microwave- 
induced Ramsey experiments on the \F = 3,mF = 0) ^ |F = 4, mi? = 0) 
clock transition of '^^Ca"'". For Ramsey interrogation periods r = 1 s, the Ram- 
sey contrast was still 88% while for r = 50 /is, a cont rast of 97% was found . 
The experimen t was carried out at a field of 0.05 mT (IBenhelm et al.l . l2008al : 



Benhelml . |2008[ ). 
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3.1.3 Imperfect control 



Other serious sources of decoherence are imperfect realizations of the intended 
gate operations, usually caused by fluctuating or insufficiently calibrated con- 
trol parameters. Typical candidates for these parameters are intensity and 
frequency fluctuations of the laser and beam pointing instabilities. In addi- 
tion, pulses on a particular ion can have unwanted side effects on the ion itself 
(AC-Stark shifts, off-resonant excitations) or on neighboring ions (addressing 
errors). Many of these errors can be greatly reduced with composite pulse 
and optimal control techniques (c.f. Sec. I2.7.3p . We list here some of these 
imperfections: 



• Pulse length errors arise from intensity fluctuations of the laser beam, 
beam pointing instabilities or just miscalibration. All of these reasons are 
almost equally relevant in current experiments. Common to all of them is 
also that the fluctuations take place at frequencies below a kHz, such that 
they can be considered constant during the execution of a pulse sequence. 
Relative amplitude fluctuations of 10~^ are typical. 

Insufficient cooling can also lead to effective intensity fluctuations. During 
each run, the phonon numbers have different values (which corresponds 
to a finite temperature of the ion crystal) and thus for a non-vanishing 
Lamb-Dicke factor the Rabi frequencies are different (see Eq. [TT]and Eq. [T3|) . 
Interestingly, for an increasing number of ions the Innsbruck group observes 
that this effect is reduced. In fact, i n their current experirn ents it is only 
relevant for single ions. As detailed by lWineland et al.l (119981 ). there are two 
reasons for this: 

(1) With an increasing ion number the Lamb-Dicke factor for each mode tends 
to get smaller. 

(2) The contributions of the increasing number of modes averages, such that 
the variance of the effective Rabi frequency narrows. 

• Detuning errors take place when the qubit transition frequency is miscal- 
ibrated or the drive frequency fluctuates slowly as compared to the duration 
of the coherent manipulation. Furthermore, magnetic fleld fluctuations have 
the same effect and usually are also slow as compared to the coherent ma- 
nipulation time. The effect of a detuning error is a constant phase evolution 
during the experiment. A simple and effective method to remove such a 
constant phase evolution is the so-called spin-echo method ( iHahru . Il950l ). 
The idea is that after half the evolution time the roles of the upper and the 
lower qubit level are exchanged. Thus the phase rewinds during the second 
half and — if the detuning is constant — arrives at zero after the complete 
evolution time. Often this method can be implemented q uite straightfor- 



wardly. An example in the ion trap context can be found in iLeibfried et al. 



( l2003bl ). However, usually during algorithms the qubits state is changed. In 
this case either more spin-e cho sequences might have to be used (e.g. for 
each free evolution one as in Barrett et al. ( 20041 )) or an effective spin echo 
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sequence can be found. For instance, iRiebe et al.l (l2004j ) optimized the time 
position of the population inverting pulses both in simulations and in the 
actual experiment. 

Addressing error While addressing a single ion with a focussed laser 
beam, residual light might affect other ions in the trap and thus perform an 
undesired unitary evolution. See Sec. 12.5.11 for more details. 
Off-resonant excitations also limit the obtainable fidelity. Off-resonant 
excitation is usually a problem if one drives a weak transition in the presence 
of a nearby strong transition (c.f. Eq. Exactl y this situati o n occ urs in 
ion traps when driving the sideband transition (jSteane et al.l . |2000| ). The 
transition matrix element of the sideband transition is weaker by a factor 
7] than the one of the carrier transition as can be inferred from Eq. [13] and 
Eq. [m Thus strong laser fields are required to obtain a reasonable gate 
speed which can then yield high gate fidelities in the presence of dephasing 
mechanisms. 

However, the strong laser field, characterized by the Rabi frequency fl, 
leads to off-resonant excitations on the carrier transition (see Eq. [9l). The 



Innsbruck experiments suffered particularly from this effect fISchmidt-Kaler et al. 



2003d ). Quantum mechanically, the off-resonant excitation can be under- 



stood as Rabi oscillations induced by a non-adiabatic switching of the en- 
ergy eigenbasis while the laser power is changed. Thus a system initially 
being in an energy eigenstate finds itself not any longer in an eigenstate of 
the Hamiltonian and consequently oscillations between the newly populated 
energy eigenstates occur. 

Off resonant excitations can be greatly reduced with pulse shapes which 
have no spectral Fourier components at the carrier-transition. In the sim- 
plest case the laser pulse powers are switched smoothly such that during 
the smooth turn on and off the system follows adiabatically. 
Unwanted AC-Stark shifts have a similar origin as off-resonant excita- 
tions (see Eq. [S]), however, affect the phase of the qubit rather than the 
population. Here the carrier transition nearby lead s to an AC-Stark shif t 
of the qubit levels and thus the qubit phase evolves (IHaffner et al.l . l2003bl ) . 
In principle, this phase evolution can be measured and then taken into 
account in the algorithm to be implemented. In the Innsbruck setup the 
problem, however, with this approach is that the acquired phase shift dur- 
ing a controlled-NOT operation is typically on the order of 207r. In order 
to obtain a reasonable phase stability of 0.1 vr, one needs an intensity re- 
producibility of better than 0.005. To achieve this intensity stability of a 
particular polarization at the ion position with a narrow beam waist is 
quite demanding. To relieve this stringent condition, a second light field 
can be used w hich induces an AC-S t ark shift of the same magnitude but of 
opposite sign (IHaffner et al.l . l2003bl : iKaplan et al.l . |2002| ). This field can be 
most conveniently derived by driving the acousto-optical modulator (AOM) 
used to control the qubit-control field with a second RF signal. Thus the 
two laser fields, for qubit manipulation and for AC-Stark shift compensa- 
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tion, are generated simultaneously by the same AOM. Both light fields pass 
then along the same path to the ions such that they pick up virtually the 
same intensity, polarization and beam pointing fluctuations, removing the 
AC-Stark shifts to a large extent. 

In the Innsbruck experiments, a lso a considerable AC -Stark shift appears 



due to dipole-allowed transitions (IHaffner et all l2003bl ). For the blue side- 
band, the presence of the dipole-allowed transitions cancels partially the 
shift induced by the presence of the carrier transitions, whereas for the red 
sideband they add. Thus in the Innsbruck experiments for coupling the ions 
to the motion, the blue sideband is preferred over the red one. Another 
possibility to reduce AC-Stark shifts is to use the polarization degree of 
freedom. The NIST group tunes the polarization of the Raman-laser beam 
pair to minimize the shift (jWineland et al.l . |2003| ). 

For Raman-driven qubits, spontaneous d ecay from the leve ls used to couple 
the two qubit levels has to be considered. lOzeri et al.l (120051 ) show that using 
very large detunings this decoherence effect can be reduced sufficiently, how- 
ever, at the expense of requiring large laser p owers. For a detailed discussion 



of these issues, we refer here to lOzeri et al.l (120071 ). 



Most of these decoherence sources can be minimized by changing external 
parameters. For example the laser intensity can be reduced such that AC- 
Stark shifts become negligible at the expense of slow gates. Slow gates in 
turn increase the susceptibility to dephasing due to fiuctuating magnetic fields 
and laser frequenci es. Similarly high trap frequencies allow for faster gates 
( ISteane et al.l . 120001 ) . but make good addressing of the individual qubits more 
difficult. In the experiments therefore often a compromise has to be made to 
keep all decoherence mechanisms reasonably small. Finally, 



The previous paragraphs listed the most common important imperfections. 
However this list is of course incomplete: for instance, a finite residual tem- 
perature of the ion string leads to single and two-qubit gate errors. Especially, 
the Cirac-ZoUer gate is very susceptible to imperfect cooling. This error source 
is very special, since the combined fidelity of two concatenated gates is not the 
product of the individual gate fidelities as discussed in Sec. 14. 7[ Therefore, a 
full simulation of the whole algorithm is required to make reliable predictions 
on its performance. 



3.2 Motional coherence 



In the Cirac-Zoller proposal, the quantum information is temporarily stored 
in the motional degree of freedom. Some other gate types, like the M0lmer- 
S0rensen gate require the ion string to be only well within the Lamb-Dicke 
regime (Sec. 12.6.31) . Both gate types are affected by population changes (mo- 
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Fig. 21. Lumped circuit model of a single ion at velocity x interacting with the trap 
electrodes. The connection between the trap electrodes is mainly characterized by 
the resistivity R and capacitance C between the electrodes. 



tional heating) and dephasing (e.g. trap frequency fluctuations) during the 
gate operation. Dephasing is often caused by slow drifts of the trap volt- 
age on the order of a few Hertz, whereas motional heating can be caused 
by electromagnetic background radiation at the trap frequencies. This back- 
ground stems predominantly from voltage fluctuations in material close to the 
trap. The most fundamental source of these voltage fluctuations should be 
the thermal motion of the electrons inside the conductors. This mechanism 
has been thoroughly iii v estiga ted, both theoretically and experimentally, by 
Wineland and Dehmeltl (119751 ) with electrons in Penning traps. The Johnson 
noise heating power P is given by: 



-Pnoise — kT/^U , 



(17) 



where kT is the thermal energy and Au is the frequency bandwidth in which 
the ion accepts the power. The time r in which one motional quantum of 
energy Eq = hv is generated is given by 



r 



-1 ^ noise 



kTAu 
hv 



kT 
hQ' 



where we introduced the quality factor Q of the ion motion. This quality 
factor can be derived from the dissipated power Pdis = I^ReZ of the current 
induced by a single ion with an energy Eion at the real part of the impedance 
KeZ. For this we assume a lumped circuit model (see Fig. [2T]) where the ion 
indu ces a current / = qx/D with q being the charge and x the velocity of the 
ion (IShockleyi . Il938l ). The characteristic dimension D is on the order of the 
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ion-electrode distance. Using this, we obtain for the quahty factor: 



Q 



Pdis/iy PReZ q^ReZ 



(19) 



Inserting Eq. [19] into Eq. [18], we arrive for the time in which one motional 
quantum is acquired at 

- ■ (2°' 

For typical values oi D = 100 /xm, ReZ = 1 Vt aX room temperature, the 
expected heating time from Johnson noise is r ~ 200 s/quantum and thus 
very small. 

In order to measure a heating rate, one can cool the ion (string) to the ground 
state, wait for some time to allow for some heating and then probe the strength 
of the motional sideband. Assuming a thermal distribution, the mean phonon 
number is deduced by employing Eq. [13] Depending on the thermal excitation, 
the Rabi oscillations on the blue sideband speed up and eventually degrade 
(c.f. Sec. 12.6.11). Repeating t his pr o cedure for var i ous w aitin g times yields the 



heatin g rate. ISeidelin et al.l (120061 ) . [Epstein et al.l (l2007l ) and lWesenberg et al. 



(I2OO7I ) developed and applied another method which is based on the strength 
of the ion fluorescence. The basic idea is that the ion motion leads to Doppler 
broadening of the absorption spectrum. The change in fluorescence is recorded 
as a function of a waiting time when cooling is switched off. The latter method 
is less sensitive than the first one, however, does not require manipulation of 
the sidebands. 



The N I ST group obsery e d hea ting rates of a few phonons per ms (iLeibfried et al. 



2003al : iTurchette et al.l . I2OOOI ) in various traps. This is orders of magnitude 



more than what is expected from fundamental electrical noise in the trap 
electrodes (Eq. [20]). Patch charges on the trap electodes have been suspected 



to ca use these excessive heating rates (IWineland et al.l . Il998l : ITurchette et al. 



2OOOI ). The former publication also discusses various other sources of heating in 
detail. Patch charges can be generated when an electron beam is used to ionize 
the atoms during trap loading. Indeed many exp eriments suggest that u sing 
photo ionization techniques to produce the ions (IKjaergaard et al.l . I2OOOI ) can 
help to reduce the heating rate. A reason for the seemingly reduced heating 
rate could be the reduced vapor pressure of the atoms during photo ioniza- 
tion as compared to the less efficient electron beam ionization. Thus a much 
reduced atom flux can be used which reduces deposition of atoms on the trap 
electrodes. Furthermore, ionization with a laser produces only a minimum of 
charged particles whereas the electron beam can charge any insulating layer 
on the trap electrodes. While photoionization seems to lead to reduced heating 
rates, no experiment with laser cooled ions was reported as of yet where the 
fundamental thermal noise given in Eq. [20] dominated the heating rate. Only 
in some Penning trap experiments using a resonance circuit to enhance the 
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Fig. 22. Single ion h eating rate as fun c tion o f ion-electrode distance (a measure of 
the trap size) (from Deslauriers et al. ( 20061 )). The trap consists of two needles to 
which the radio frequency is applied to. Furthermore, two heating rates from planar 
surface traps are added: the cross shows the heating rate of a ^^ Mg"*" ion in a room 
temperature trap (electrode material: gold) ( Epstein et al. . 200?! ). while the triangle 
shows t he heating rate of a single Sr"*" ion where the silver trap electrodes were held 
at 6 K (jLabaziewicz et al.l . l2008al ) . 



real part of the resistivi t y in Eg. [201 thermal noise plays the dominant role 
(IWineland and Dehmeltl . 1 197.4 iHaffner et all . bnO.Sal ). 



The Ann-Arbor group o bserved much reduced h eating rates by cooling the 
trap electrodes (see Ref. iDeslauriers et al.l (120061) and Fig. [22]). In the course 
of cooling the trap electrodes from 300 K to 150 K, the heating rate dropped 
by more than one order of magnitude. This strong dependence on the temper- 
ature hints at a thermally activated process causing the unexplained heating 
in ion traps. In addition, the authors measured heating rates as a function 
of the trap size d. From Eq. [201 ^ l/cP scaling is expected while one ex- 
pects a dependerice in the case of hea ting due to moving patch charges 
(iTurchette et all . I2OOOI : lEpstein et aLl . l2007l ). T he Ann- Arbor group extr acted 
from the data in Fig. [22] an exponent of 3.5 (iDeslauriers et al.l. 120061) on ly 



slightly different from the postulated exponent of 4 (iTurchette et al.l . I2OOOI ) 



The MIT-group has investigated various planar traps made of silver electrodes 
on a quartz substrate clos e to 4 K as well as one trap at room temperature 
(jLabaziewicz et al.l . l2008al ). They measured a heating rate as low as 2 quanta/s 
for trap sizes on the order of 100 /zm while a similar trap at room temperature 
had a devastating heating rate seven orders of magnitude larger. Even the 
extremely small heating rates for the 4 K experiments cannot be explained by 
Johnson noise, only. They also measured heating rates for three different trap 
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sizes and found them to be consistent both with a l/rf^ and a l/rf^ scahng. In 
addition, a strong dependence of the heating rate on the anneahng tempera- 
ture used in the fabrication process was found. Furthermore, the MIT group 
investiga te the heating rate of the i ons as a function of the electrode temper- 



ature T (ILabaziewicz et al.l . l2008bl ). Above T = 40 K, they found that the 



heating rate is proportional T^^ with 2 < /3 < 4 depending on the trap. Over- 
all, these findings suggest that with improved fabrication methods, smaller 
heating rates can be achieved. Furthermore, the NIST group observed heating 
rates of 300 quanta/s for a ^^Mg"*" ion 40 /im above the gold surface of a planar 
trap (lEpstein et al.l . 120071 ). This heating rate is significantly smaller than what 
one would expect from the MIT measurements and supports the conclusion 
that choice of materials and fabrication methods are very important to achieve 
small heating rates. 



To speed up quantum gates and to ease cooling, there is a strong tendency 
towards constructing small ion traps. On one hand, small traps with charac- 
teristic dimensions of a few tens of microns allow for large trap frequencies 
on the order of 10 MHz, on the other hand they seem to lead to inacceptable 
heating rates. Therefore, heating in ion traps is not only an interesting topic 
on its own but needs a lot of attention from a technological point of view. 
Cooling the traps to 4 K seems to offer a solution to the heating problem. 

Many two-qubit gate implementations, however, store quantum information 
in the motional degree of freedom. Therefore, dephasing of the motional states 
must be also taken into account. To measure the motional coherence, a su- 
perposition of two motional states can be created whose phase coherence is 
tested after some waiting time. In the Innsbruck experiments, the pulse se- 
quence (read from right to left) R'^^ti, 0)i?'-^(7r/2, 0) creates the state 0) + 
\D,1))/V2. The inverse pulse sequence R'-'{n/2,ip)R'^{7r,0) closes the inter- 
ferometric procedure after some waiting time T. From the contrast of the 
interference fringes obtained by varying (f with waiting time T, the coherence 
time can be directly deduced. One might be tempted to use just a pair of two 
R^{tt/2) pulses to implement the Ramsey experiment. In this case, however, 
motional as well as electronic dephasing mechanisms lead to decoherence of 
the intermediate state (l^, 0) + \D, l))/\/2. The former pulse sequence, how- 
ever, is insensitive to phase decoherence of the electronic qubit, and waiting 
times T of many tens of milliseconds are possible. Thus, with the former pulse 
sequence the trap frequency can be easily determined with an accuracy of a 
few Hertz. 



Using this method, th e Innsbruck group ob s erved a motional coherence time 
on the order of 100 ms (jSchmidt-Kaler et al.l . l2003al ) on the center-mass-mode, 
being consistent with expected voltage fluctuations on the order of 10~^. Fur- 
thermore, the Oxford group observed on a single ^°Ca+ion a motional coher- 
ence time of 182 (36) ms, limited most likely by motional heating of about 
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axial 



Fig. 23. An anti-correlated radial motion (excited rocking mode) of a two-ion crystal 
changes the mean distance between the two ions. 

3 quanta/s ( Lucas et al.l . 20071 ). However, the Innsbruck group found that for 
the axial breathing mode (a nd other hig her axial modes) coherence times of 
about 5 ms are more typical (lRoosl . l2008l ). Thus it must be concluded, that on 
the breathing mode, a depha sing raechaii ism is present which can be neglected 
for the center-of-mass mode. iRood (120081 ) show, both theoretically and exper- 
imentally, that for a two-ion crystal the breathing mode frequency depends 
on the motional state of some of the radial modes. The basic mechanism re- 
sponsible is that an anti-correlated motion (the rocking mode) along a radial 
direction changes the mean distance between the ions (see Fig. 1251) . Thus, the 
repelling force between the ions is reduced leading to a reduction of the ax- 
ial breathing mode frequency with increasing excitation of the rocking mode. 
Therefore, if the rocking modes are not in a well-defined state, the breathing 
mode frequency is different for each experimental realization which in turn is 
interpreted as dephasing of the breathing mode. 



3.3 Modelling ion trap quantum computers 



The influence of each imperfection on the performance of the quantum com- 
puter can be reliably estimated with numerical simulations. In this way the 
necessary steps to improve the performance can be analyzed and planned in 
detail. For ion trap quantum computers this procedure i s relatively straightfor- 



ward as the Hamiltonians are well-known (se e Eq. [2]) (IWineland et al.l . Il998 



Leibfried et al. , 2003b : Jonathan et al. , 2000l ). Here we describe the simula- 
tions as appropriate for the Innsbruck set-up. With small modifications these 
simulations should also be applicable to other ion trap set-ups. 

Eq. [2] is a good starting point to model ion trap quantum computers. The 
Hamiltonian is first generalized to multiple ions: 



(21) 



Here the indices n and m denote the various ions and motional modes taken 
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into account, respectively, and rinm accounts for th e different coupling strengths 
of the ions to the motional modes (jjamesl . Il998l ). 



For the simulations the initial state vector is first transformed into the rotat- 
ing frame of the laser field such that Eq. [21] becomes time independent and 
the Hamiltonian can be directly integrated. A quantum algorithm usually con- 
sists of laser pulses with different frequencies. Therefore, this procedure has 
to be carried out for each pulse, separately. For the bi-chromatic laser fields 
employed for the AC-Stark shift compensation (c.f. Sec. 13.1.31) and M0lmer- 
S0rensen-gates (Sec. 12. 6p this method fails. In these cases the differential equa- 
tion can be numerically integrated. In the following we describe how almost 
all experimental imperfections were incorporated into the simulations: 

AC-Stark effects due to the carrier transition are described by Eq. [51] and 
appear with increasing fi„ (see Sec. 12. 2p . Therefore, these shifts are auto- 
matically included in the simulation. AC-Stark shifts due to other (dipole) 
transitions are not taken into account by Eq. [21] However, in the experi- 
ments a second light field is used anyways to minimize the total effect of 
all AC-Stark shifts. Therefore the light shifts must be artificially removed 
from the Hamiltonian to match the experiments. 

Off-resonant excitations (c.f. Eq. [9]) are described by Eq. [2T]and are there- 
fore automatically included. 

Laser freqency noise and magnetic field noise are divided into fluctu- 
ations slower and faster as compared with typical coherent manipulation 
times (on the order of 1 ms in current experiments). In the experiments 
usually slow fluctuations dominate and the laser detuning remains to a large 
extent constant during coherent state manipulation. The observed increase 
in the coherence time — when spin echo sequences are used — supports this 
conjecture. These slow fluctuations can be modeled by running the simula- 
tions for several detunings from the qubit transition and by averaging the 
measured populations. On the other hand, fast fluctuations can be taken 
into account by transforming Eq. [21] into a master equation. In this case, the 
dimension of the state space describing the system is squared as compared 
to the Schrodinger approach and thus this method starts to get tedious 
already for simulations of a five-ion algorithm. In fact, even today's super- 
computers cannot hold the complete density matrix of an arbitrary twenty 
qubit system in their memory. 

Laser intensity fluctuations are assumed to be constant during the coher- 
ent state manipulation and therefore can be modeled by simulating the 
algorithms a few times for various laser intensities. 

Addressing errors are described by setting the Rabi frequencies f2„ in Eq. [51] 
to the corresponding values. However, it should be noted that there is an 
additional degree of freedom connected to addressing errors: the phase of 
the laser at an ion position when directed on this particular ion and when 
mainly directed on another ion is in general not the same. This is due to 
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the different paths the hght field takes in these two cases. Usually, in the 
simulations all phases are taken to be the same as they seem to give an 
upper bound for most algorithms. 
Imperfect ground state cooling can be taken into account by running the 
simulation with the different initial states and averaging the results appro- 
priately. 

These simulations allow one to reduce the sensitivity of quantum algorithms to 
the respective imperfections. In the Innsbruck experiments, especially the in- 
fluence of laser freqency and addressing errors issues was reduced considerably 
by optimizing the implementations taking the simulations as a guide. 

Finally, we document some of the imperfections not taken into account in the 
simulations for the Innsbruck experiments. In Eq. [21] each ion was approxi- 
mated as a two level system. For ^°Ca"^ the 5'i/2-ground state is split into two 
Zeeman levels and the upper qubit level £^5/2 into six levels. The level sepa- 
ration is about 5 MHz and therefore not much larger than the involved trap 
frequencies. Therefore, off-resonant excitations of the additional transitions 
are possible. Furthermore, AC-Stark effects arise due to the presence of th e 



other Zeeman levels as well as due to dipole transitions (IHaffner et al.l . l2003bl ). 
Finally, the AC-Stark effect was canceled with a second off-resonant laser field. 
Treating this second light field leads to a time dependent Hamiltonian and was 
therefore usually not taken into account in the simulations. 



4 Key experiments 



4-1 Cirac-Zoller-type gates 



The NIST group demonstrated the central part of the Cirac-Zolle r gate on a 



single ^Be"*" ion by implementing the operations displayed in Eq. [T^ (IMonroe et al. 



1995al ). The phase gate was turned into a controUed-NOT operation by insert- 
ing two Ramsey 7r/2 pulses, one before and one after the phase gate, to verify 
the quantum nature of the phase gate. The first two motional excitations 
n = {0, 1} (see Fig. [H ciJtrap = 27r x 11 MHz) served as the control bit, while 
the target bit was represented by superpositions of the \F = 2,mf = —2) and 
\F = l,mf = —1) states. For the auxiliary state, the \F = 2,mf = 0) state 
was used. 



The co mplete Cirac-Zoller gate was finally implemented by lSchmidt-Kaler et al. 



( l2003d ): two ^^Ca"*" ions were addressed individually using a laser beam with 
a waist of about 2/im. First, the quantum state of ion 7^1 was mapped onto 
the breathing mode with a blue sideband pulse i?"'"(7r). Then, a controUed- 
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NOT gate with the motion as the control bit and the electronic state of 
ion #2 as the target bit was carried out, before the motional state was 
mapped back onto ion ^1. However, to implement the conditional phase 
shift, the Innsbruck group did not use an auxiliary level but the compos- 
ite pulse sequence R^(Til2,Ti I2)R^{tt I ^/2,Q) R+(ti /2,ti I2)R+(ti I ^/2 , ti) dis - 



cussed in Sec I2.7.3I In the first experiments (jSchmidt-Kaler et al.l . l2003d ) 



the controlled-NOT mapped the four logical eigenstates and one superposi- 
tion state to their desired respective output states with an average fidelity 
of 0.73(2) (coherence time r ~800 /xs, addressing error e ~0.05). After sev- 
eral improvements of the experimental set-up (r ~2 ms, e ~ 0.03, more 
f lexible computer c ontrol, etc.), fidelities as large as 0.91.0(6) were observed 



ipi 

( iRiebe et al.l . l2006l ) where here the fidelity is defined as an average of the 



overlap of the produced output states with the ideal output. Additionally, 
this work implemented the conditional phase shift with the pulse sequence 
i?+(7r/2, 0)i?+(y27r, 7r/2)i?+(7r/2, tt) yielding a fidelity of up to 0.926(6). 



4.2 Entangled states with trapped ions 



One important application of quantum computers in basic research is the gen- 
eration of interesting quantum states as for insta nce the first determiri istic 



generation of entangled particles with ^Be"*" ions (ITurchette et al.l . Il998l ). In 



addition, ne w applications of entangled states especially for metrology appear 
constantly (IBlatt and Winelandl . l2008l ). Entangled states play an important 
role in discussions on the foundations of quantum mechanics. Especially since 
Bell formulated inequalities which coul d dis t inguish between so-called local re - 



alistic theories and quantum theories (IBelll . Il965l . Il97ll : IClauser et al.l . Il969l ) 



physicists were keen to produce these states and to check the predictions 
of quantum mechanics. Since then, there have been numerous experiments 
demonstrating a violation of a Bell inequality. Almost a ll of these experiments 
were carried out wi t h pho tons (for a sumrnary se e e.g. IClauser and Shimony 
fll978h . lWeihs et al.l (Il998h and lTittel et~aD (Il998h l Photons naturally explore 
the non-locality of entanglement and thus violations over distances of many 
kilometers were established. However, current detection efficiencies of photons 
are not high enough to close the so-called 'detection loophole', i.e. one must as- 
sume that the detected particles rep resent a fair sample of all particles emitted 



by the source (jClauser et al.l . Il969l ) . Therefore there is a large interest in test- 



ing Bell- inequalities with trapped ions where the detection fidelities approach 
unity. Such an experiment was conducted with ^Be+ ions by lRowe et al.l (|200ll ) 
and closed the detection loophole. However, the ions were not detected outside 
their respective light cone, i.e. the detection time was longer than the time it 
takes light to travel between the ions. Thus there could still exist a combined 
detection-locality loophole. This loophole could be excluded for instance by 
creating entanglement between ions or atoms separated by several kilometers. 
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Another interesting aspect of the ion trap experiments is that the entangled 
states are produced deterministically. That means that in contrast to the pho- 
ton experiments, the entangled states are created on demand. More impor- 
tantly, the entangled states are not destroyed during their creation such that 
they can be used for further experiments. 

Current research is directed towards entangled states with more than two 
particl es. Already for three qub its, two classe s of entang l ement appear: GHZ 



state s (iGreenberger et al.l . ll989l ) and W states (iDiir et al.l . l2000l : IZeilinger et al. 
I992I ). These two classes of entanglement cannot be transformed into each 



other by local operations and classical communicatio n, i.e. with single -qubit 



operations and measurements of the individual qubits (iDiir et al.l . I2OOOI ) . Both 



classes are not only maximallay entangled but also violate Bell-type inequali- 
ties. 

GHZ states are states of the form 



I GHZ, 



I00---0) + |11---1))/V2. 



f22) 



GHZ states with many qubits can be interpreted as Schrodinger-cat states. For 
this, e.g. the first qubit is treated as a separate degree freedom and all other 
qubits as a single quantum system, i.e. the cat system. Then the state of the 
first qubit indicates the state of the second "macroscopic" system. Another 
use of GHZ states is by encoding quantum information in a superposition of 
the form a|000) -|- If a single qubit flips, the original st ate can still 



be recovered with so-called quantum error correction protocols (jShorl . 11995 
Steanel . [T996h . 



W states are states of the form 



|Wn) = (|0 ■ ■ ■ 001) + |0 ■ ■ ■ 010) + |0 ■ ■ ■ 0100) + ■ ■ ■ + |10 ■ ■ ■ 0)) /Tn . (23) 



They are remarkably stable against various decoherence s ources: they are in - 
trinsically stable against collective dephasing mechanisms ( Roos et al. . 2004bl ) 
and even loss of qubits does not completely destroy the entanglement present 
in them. 



A four ion GHZ state was first produced by iSackett et al. fl2000f ) usmg a 
M0lmer-S0rensen type gate (see Sec. I2.6.3p . The produced GHZ-state was 
analyzed by applying i?'-^(7r/2, y9)-pulses to all ions simultaneously and mea- 
suring the number of fluorescin g ions with a photo multiplier. Here, we will 



illustrate this procedure used by ISackett et al.l (l2000l ) for the two particle Bell 
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state (|00) + |11))/V2: 

Rf(7r/2,(^),Rf(7r/2,</.) 



|00) + 111) 



(|0) +ie^'^|l)) (|0) +ze*^|l)) + (|1) +ie-^'^|0)) (|1) +ze-^'^|0)) 
= (1 - e-2*^)|00) + ie'^il + e-2^^)|01) 
+ ie'^{l + e-2^'^)|10) + (1 - e-2^^)|ll) , 

(24) 

where we used Eq. (TUlin the first step and omitted the normalization factors. 
To find an estimate for the fidehty of the original Bell state, it is useful to 
introduce the parity operator P which is defined as P = Pqo ~ -Poi ~ -Pio + -Pii- 
Here P^^y are the probabilities to find the ions in state \xy). Plotting the 
expectation value of the parity, we see that it oscillates twice as fast as the 
phase if of the analyzing pulses. This behavior can also be interpreted as a 
consequence of the doubled energy difference between the |00) and the |11) 
state as compared to the single ion case. 

To find the fidelity of the Bell state, two sets of experiments can be carried 
out: first a Bell state is created and the populations Pqo and Pn are recorded. 
In a second set of experiments, the maximum (max P) and minimum (min P) 
of parity oscillations as described in Eq. [2l] are determined. The overlap F of 
the experimentally produced state with a state of the for m (|00) +e'''''|ll))/\/ 2 



is given then by P = (Poo + Pii)/2 + (maxP-minP)/4 fISackett et all . 120001 ). 
No individual addressing of the ions is required in the analysis procedure. In 
addition, the parity can be inferred from the global fluorescence of the ion 
string. Thus, this analysis method is relatively simple and efficient. Further- 
more, it can be generalized to GHZ states with an arbitrary number of ions and 
is thus very useful to gain information on the generated GHZ states, without 
individual qubit addressing and read-out. 

This analysis technique was also used to verify the creation of a three particle 



GHZ-states with fidelities of up to 0.89 (ILeibfried et al.l . 120041 ). The quality of 
the GHZ state was high enough, that from the resulting generalized Ramsey 
fringes (c.f. Eq. [2^ the phase could be estimated 1.45 times more accurately 
than using three uncorrelated particles. Estimating the phase of superposi- 
tions is quite important in frequency measurements. The high fidelity of the 
GHZ states was made possible by using the geometric phase gate (Sec. I2.6.4p . 
Encouraged by this, the NIST group applied this gate also to create four-, 
five- and six-particle entangled GHZ-states, with lowe r bounds for th e fidel i- 



ties of 0.76 (1),0.60 (2) and 0.509 (4), respectively (ILeibfried et all . [20051). 
As for GHZ states, a fidelity above 0.5 implies automatically the presence of 
genuine A^-partite entanglement, the latter experiments demonstrated up to 
six partite GHZ-like entanglement. 
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Fig. 24. Pulse sequence to generate a GHZ state. The controlled-NOT operation 
is implemented with a composite phase gate $(7r) (see Sec. I2.7.3P sandwiched in 
between two i?'-'-pulses on the center ion. Hatched areas indicate sideband pulses. 

While the NIST group uses predominantly global addressing and state read- 
out, the Innsbruck group entangled ions mainly with laser pulses addressed 
to individual ions. For instance, a Bell can be creat ed in the following way 
(please note that the right-most ion is the first one) (IRoos et al.l . l2004al ): 



"^'-^^ {\DS,Q) + e'^\DD,l))/V2 



(25) 



The success of the Bell-state generation is usually verified using a procedure 
called state tomography (see Sec. 14.41) . The laser phase offset ip of the first pulse 
determines phase of the Bell state. In addition, an additional i?^(7r, 0)-pulse 
on the second ion transfers the (l^^, 0) + e''P\SD, state to {\SS, 0) + 

e^^|DD,0))/V2. Thus, using this toolbox all four Bell-states can be created 
in the same set-up. 

Furthermore, the Innsbruck group used the fl exibility of the entan gling method 
to create three partic le GHZ and W s t ates (IRoos et al.l . l2004bl ). The idea to 
create a GHZ state dGirac and ZoUej . Il995l ). is to apply a controlled-NOT 
gate while the motional degree is still in a superposition of |0) and |1) (see 
Eq. [25]) . For this CNOT-operation, the motion i s the control bit and the new 
ion is the targe t bit (t he second line in Fig. [24]) (jRauschenbeutel et al.l . l2000l : 
Raimond et al.l . l200ll )). Inserting more and more CNOT's, this GHZ-state 
generation method is straightforwardly generalized to more particles. 

In order to create a three-ion W state IW3) = (|DDS + |DSD) + |SDD)))/y3, 
Eq. [25] can be generalized differently: The length of the first blue sideband 
pulse is adjusted such that the state (jS'S'S', 0) + \/2\SSD, \))/\/?) is created 
and ions #2 and #3 are flipped to obtain {\DDS,Q) + y/2\DDD,l))/^/?,. 
Then the phonon is shared between the remaining ions #2 and #3 to create 
{\DDS, 0) + \DSD, 0) + \SDD, 0))/V3. Fig[25] shows the corresponding pulse 
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ii+(2arccos(l/\/3),0) 

,///// / 









i?+(7r/2,7r> 



R^'iTT, 0) 



Fig. 25. Pulse sequence to create the W state i\DDS) + \DSD) + \SDD))/^/3 from 
the \SSS) state. Hatched areas indicate sideband pulses. 



sequence. It can also be directly generalized to ions by adjusting the pulse 
area of the first beam splitter to arccos ( 1/\/N) and th e n shar ing the phonon 
excitation among the other ions equally (IHaffner et al.l . l2005al ). Interestingly, 
the required blue sideband pulse area grows only logarithmically with the ion 
number and thus the generation time grows sublinear with the number of 
ions. Therefore, this scheme opens the possibility to generate large entangled 
states. I n experiments, W stat es ranging from four up to eight ions have been 



created (IHaffner et al.l . l2005al ). Using a technique called state tomography 



(see Sec. 14.41) . the experimentally obtained states have been fully charactized 
(Fig. HTj). Analysing the measured density matrix showed that the generated 
states indeed carried genuine iV-particle entanglement. 



We add that an even faster scheme to produce W states was pro posed which 



does not require ind ividual addressing in the entangling procedure (IRetzker et al. 



20071 : ISolanol . l2005h . The idea here is that first a \DD ■ ■ ■ D,l) state is created. 
Then a R^^ti, 0) pulse addressed to all ions is supposed to generate the desired 
W state by mapping the phonon to the electronic state of one of the ions, i.e. 
creating a symmetri c superposition with exactly one electronic state flipped. 
Retzker et al.l (120071 ) further generalized this procedure to W states with more 
than one excitation (Dicke states). 



Entangled states h ave also been produce d with trapped ^^^Cd^ a nd ^^Ca"*" ions 
by th e Ann- Arbor ( Hal] an et al. . 2005bl ) and the Oxford groups ( Home et al. . 
20061 ). respectively. For the two qubit levels, the Ann- Arbor group used the 
\F = 0,mp = 0) and \F = l,mp = 0) state of the ground state of ^^^Cd"*", 
taking advantage of its insensitivity to the Zeeman effect in first order. The 
geometric phase gate ( Sec. I2.6.4p doe s not work efficiently on magnetic field 
insensitive transitions (iLangerl . l2006l ). Instead, the Ann- Arbor group used a 



M0lme r- S0rensen gate (Sec. 12.6.31) to entangle the two ions (iHaljan et al. 



2005al lbh. Furthermore, a tomographic state characterization was applied to 
evaluate the produced states and the degree of entanglement thoroughly (see 
Sec. 14. 4p . To achieve the required individual addressing capability, a combi- 
nation of ion selective AC-Stark shifts and microwave fields was used (c.f. 
Sec. l^XTj) . 
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The O xford group created Bell states with two trapped ^^Ca"*" ions (iHome et al. . 



20061 ) ■ They encoded the quantum information in the Zeeman manifold of the 
S'i/2-ground state, thus effectively using the direction of the valence electron's 
spin. Adopting the geometric phase gate (Sec. I2.6.4p for ^^Ca"*", they created 
the Bell state (|TT) - |ii))/v^. 

In Innsbruck, a M0lmer-S0rensen gate was used to prepare a pair of ^^Ca"*" ions 
in the entangled state if) = {\SS) + i\DD))/\/2, where \S) = \Si/2,m = 1/2) 
and \D) = \D5/2,m = 3/2). Using a gate time r = 50 /is with the laser 
light being smoothly swit ched on and of f within 2.5 /is, a Bell state fidelity 



of 0.993(1) was achieved (IBenhelm et al.l . l2008bl ) when the ions were cooled 



to the ground state of the motional mode mediating the coupling. Moreover, 
uneven multiples A; = 1, 3, 5, . . . , 21 of the gate were used to create entangled 
states. For k = 21, the state fidelity was still 0.8. Finally, the gate yielded 
Bell states with a fidelity of 0.96 even with the ion string cooled only to the 
Doppler limit (< nbus >= 17). 



4-3 Decoherence free subspaces 



Laser frequency and magnetic field fluctuations are usually the dominant de- 
coherence mechanisms in ion traps. Both mechanisms lead to fluctuations of 
the phase between the laser and the atomic polarization and thus to dephasing 
of each qubit, however, to a good approximation by the same amount for all 
qubits. If one encodes a single qubit in two ions in such a way that the two 
phase evolutions cancel each other, the original qubit is protected from this 
global dephasing and the quantum information is encoded into a decoherence 
free subspace (DFS). In particular, superpositions of the form a|01) + /3|10) 
are transformed by the global single-qubit phase-change |1) —>■ e*'^|l) in the 
following way: 

a|01) + /5|10) ^ ae^'^lOl) + (3e''''\10) . (26) 

The global phase factor e*"^ cannot be observed, such that the state remains 
immune against collective dephasing. 



This property was demonstrated by iKielpinski et al.l (120011 ) using an engi- 
neered dephasing mechanism. The qubits were e ncoded in the hyperfine -states 
of ^Be"*". As a controlled dephasing mechanism, IKielpinski et al.l (120011 ) chose 
an unfocused off- resonant laser beam with "random" intensity. The laser beam 
leads for each experimental realization to a different AC-Stark effect, however, 
common to both ions. 



In an environment with natural dephasing, iRoos et al.l (l2004al ) observed also 
much increased lifetimes of qubits encoded in a DFS as compared to single ion 
qubits. In these experiments, the coherence of a state formed by two ^"Ca"*" ions 
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Fig. 26. Energy level diagram for two ions with ground state |0) and excited state 
A collective-dephasing free qubit is formed by the degenerate logical basis 
{|10),|01)}. 

I^qubit) = {\DS)+e^'-^\SD))/^/2wa,s retained for 1 s, the dominant decoherence 
mechanism being spontaneous decay of the -D5/2 level. For comparison, laser 
frequency noise and magnetic field fluctuations led to a single ion coherence 
time of 1 ms. 



Furth ermore, both the Innsbruck f Haffner et al. . 2005bl ) and the NIST ( Langer et al. 



20051 ) groups encoded quantum information in the ground states of two-ion 
strings with ^^Ca"*" and ^Be"*", respectively. Coherence times of 34 s and 7 s, 
respectively, were measured. In both cases, fluctuations of the magnetic field 
gradient were believed to be the reason for the decoherence in the DFS. For 
these experiments, extreme care must be taken to switch off the laser light 
properly. Not only residual light scattering rates on the order of 0.1 photon/s 
destroy the coherence, but also fluctuating differential AC-Stark shifts on the 
order of 1 Hz destroy the phase coherence of the entangled states. Finally, we 
note that these experiments demonstrated also extremely long lived entangle- 
ment of up to 20 seconds between two parties separated by 5 fim. 

Apart from robustly encoding quantum information, decoheren ce-free sub 



spaces have also found an application in quantum metrology. In iRoos et al. 



( I2OO6I ). a Bell state was encoded in a combination of Zeeman sub levels of the 
-D5/2 level of two ''^Ca^ ions. The state was decoherence-free with respect to 
fluctuations of the magnetic field but sensitive to energy level shifts caused 
by static electric field gradients. In this way, the quadrupole moment of the 
metastable state could be determined with high precision by monitoring the 
Bell state's phase evolution over a duration orders of magnitude longer than 
the single-qubit coherence time. 



4.4 State tomography 



Quantum state tomography (jParis and Rehacekl . l2004l ) is a measurement tech- 
nique that provides access to all the information stored in density matrices 
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describing pure and mixed quantum states. It requires the quantum state of 
interest to be available in rnany co pies. While the basic measurement principle 
dates back fifty years (IFanol . 119571 ). experim ental implemen t ations of quanturn 



state tomography started o nly in the 1990 's (jSmithey et al.l . Il993l : iDunn et al. 



19951 : iLeibfried et al.l . 119961 ). Tomographic measurements of systems composed 
of qubits have been implemented in experiments with nuclear magn e tic res- 



onance, pho t ons, t rapped ion s and superconducto r s (IChuang et al.l . Il998a 



White et al.l . Il999l : iRoos et al.l . l2004al : ISteffen et all . 120061) 



Noting that the density matrix of a single qubit can be represented by 



(27) 



we see that the density matrix of a single qubit can be inferred by measuring 
the expectation values (o"q,), (a = x,y,z), of the Pauli spin matrices. The 
measurement of az is accomplished by projecting the qubit onto its energy 
eigenstate basis. For the measurement of (Tx,y, an additional 7r/2-pulse of suit- 
able phase precedes the projective measurement. The tomographic procedure 
can be easily extended to systems of several qubits by measuring the joint spin 



expectation values a, 



(ni), 



("2) , 
02 



0-^k) where an 



denotes a spin component 



of qubit nj{aaj € {I ,ax,(Jy,(Jz})- This way, the determination of the density 
matrix of an N-qubit system requires the measurement of 4^ expectation val- 
ues. As some of the operators commute, a total of 3^ measurement bases is 
necessary. While in principle the number of measurements could be re duced 
by projecting onto mutually unbiased bases (IWootters and Fieldsl . 119891 ) . this 
procedure is of no practical importance in current ion trap experiments as 
it would demand high-fidelity entangling gate operations for mapping the re- 
quired bases to product state bases. 

A slight complication arises since in every experimental implementation of 
quantum state tomography, expectation values are never exactly determined 
but only estimated based on a finite number of measurements. The naive 
replacement of the expectation values (cTq,) in Eq. [27] can give rise to unphys- 
ical density matrices with negative eigenvalues. This problem is avo ided by 
employing a maxi r num likelihood estimation of the density matrix (IHradill . 
I997I : IJames et al.l . I2OOII ) that makes use of the estimated expectation values 
for searching in the set of meaningful density matrices the 'most likely one'. 
The maximum likelihood algorithm identifies the sought-for density matrix 
with the one that maximizes the probability of observing the experimentally 
recorded set of measureme nt results. Everi though maximum likelihood esti- 
mation has been criticized (iBlunie-Kohoutl. 120061) for bein g less accurate than 
Bayesian estimation techniques (jParis and Rehacekl . l2004l ). it has the practical 
merit of being easily implemented in experiments. 



Starting with IRoos et al.l (l2004al ). almost all Innsbruck experiments made 
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Fig. 27. Abs o lute va lues of the density matrix of an eight-ion W state (from 
Haffner et all (|2005al l). In addition, in the upper right corner the CCD-image of 
an eight-ion string is displayed. For the state tomography, the experiment was re- 
peated 100 times for each of the 6561 measurement settings. The total measurement 
time amounted to more than 10 h. The fidelity of the W state was determined to 
0.722(1). 

heavy use of quantum st ate tomography. Even an eight particle W state has 



been fully characterized (IHaffner et al.l . l2005al ). Figure [271 shows the experi 



mentally obtained density matrix of an eight-ion W state. A particular merit 
of quantum state tomography is that all physically available information on 
the quantum register is extracted. Thus all aspects of the generated states can 
be thoroughly analyzed without taking new data. 

An alterna tive to the above d escribed method, was demonstrated by the Ox- 
ford group (iHome et al.U2006l ). They used a refined, albeit partial tomographic 
procedure: instead of choosing three measurement settings for each qubit (ei- 
ther one of the pulses i?^(0, 0), i?^(7r/2, 0), or i?^(7r/2, 7r/2) preceding the 
qubit detection), they choose to apply before the detection the pulses {6i, ip) 
with 9i either 0, O.GGtt, or 0.547r and evenly distributed phases G [0,27r[. In 
this way the reconstruction is less biased and thus more robust against sys- 
tematic errors or equivalently against decoherence. 



4-5 Selective read-out of a quantum register 



For some quantum algorithms like teleportation and most error-correction 
protocols a part of the quantum register has to be read out while leaving the 
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Ion #3 




Fig. 28. Partial level scheme of the three Ca-ions ffrom IRoos et al.1 (|2004bl l1. Only 
ion ^3 is read out. Ion #1 and #2's quantum information is protected in the Zeeman 
manifold of the D5/2-level, namely the mj = —1/2 and mj = —5/2 levels. 



rest of the re gister intact. Both th e Innsbruck and the NIST group succeeded 
in this task. (IBarrett et al.l . |2004| ) employed segmented traps to separate the 
ions to be read out from the ions which should remain coherent. Now one set 
of ions can be illuminated safely with detection light while the other ions are 
left dark. 



The Innsbruc k group chose a diff e rent route to selec tively read out the quan- 



tum register (IRoos et al.l . l2004bl : iRiebe et al.l . 120041 ). Qubits were protected 



from being measured by transferring their quantum information to superposi- 
tions of levels which are not affected by the detection, that is, a light scattering 
process on the S1/2 Pi/2-transition in Ca+. In the experiments, a vr pulse on 



the S 



1/2 



•'5/2 {t^j = — 5/2)-transition transfers the quantum information 
into the {D' = D5/2 (m^- = -5/2), D = D5/2 (m^ = -1/2)} manifold. Fig. [28] 
shows two ions which are protected from the detection light at 397 nm and 
the third ion with the original encoding which is measured. After the selec- 
tive readout, a second set of vr-pulses on the D' to S transition transfers the 
quantum information back to the original computational subspace {D,S}. 



It is interesting to apply the selective read-out to an entangled qubit register 
and to d emonstrate the coll apse and even partial collaps of a wave function. 
For this IRoos et al.l (l2004bl ) first prepared a three-ion GHZ- and a W-state 
and then detected one of the ions while the quantum information of the other 
ions was still protected in the D-level. Fig. [22] shows the results of these mea- 
surements. The quantum nature of the GHZ-state was completely destroyed 
by measuring a single constituent, i.e. it was projected into a mixture of |5'5'5') 
and \DDD) (Fig. [29b and Fig. [29b). By contrast, for the W-state, the quan- 
tum register remained partially entangled as coherences between ion #1 and 
#2 persisted after measuring ion #3 (Fig. [2% and Fig. 12^). 



62 





Fig. 29 . Absolute values of density matrices after measuring ion ^^3 (from lRoos et alj 
( 2004bl )). (a) shows the density matrix of a GHZ-state before measuring and (c) 
after ion ^^3 is measured. The same for a W-state ((b) before and (d) after the 
measurement of ion #3). 



4-.6 Conditional single-qubit operations 



One can take the partial read-out of a quantum register one step further and 
perform operations conditioned on the read-out result. As will be discussed 
in Sec. 15. 2[ both the NIST and the Innsbruck group der nonstrated this pro- 



cedure with i n thei r respective teleportation experiments (iBarrett et al.l . 12004 



Riebe et al.l . |200J). Furthermore, the Innsbruck group employed conditional 



operations to deterministically transfer a three-parti c le GHZ -state with local 
operations into a two-particle Bell state (iRoos et al.l . l2004bl ). This procedure 
can also be regarded as an implementa tion of a three-spin quantum eraser as 
proposed by lGaristo and Hardyl (119991 ). 



In the experiment by iRoos et al.l fl2004bf ). first the GHZ-state i\DSD) + 
\SDS))/\/2 was created. Application of Rs{Tc/2,3n/2) yielded the GHZ state 
\D){\SD) - \DS)) + |^)(|5D) + \DS))/2. Measuring ion #3, projected ions 
#1 and #2 either onto {\SD) - \DS))/V2 or onto {\SD) + \DS))/y/2 with 
ion #3 indicating in which of the two states the first two ions were (see 
Fig. [5Uk). This mixture of the two Bell states can then be transferred to 
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Fig. 30. Density matricies during the individual steps of the determ inistic gen- 
eration of a BeU state from a GHZ-state (from Roos et al. ( 2004bl )). (a) Real 
part of the density matrix elements of the system after ion 7^1 of the GHZ-s- 
tate {\DSD) + \SDS))/V2 has been measured in a rotated basis, (b) Transfor- 
mation of the GHZ-state {\DSD) + \SDS))/\/^ into the bipartite entangled state 
IS") (IDS') + \SD))/y/2 by conditional local operations. Note the different vertical 
scaling of (a) and (b). 

a pure Bell state by inducing a phase shift of tt on ion #2 (pulse sequence 
i?^(7r, 7r/2)i?^(7r, 0)) if, and only if, ion #3 was measured to be in the D-state. 
In addition, the state of ion #3 was reset to IS*). Figure 1301 shows the interme- 
diate result before applying the conditional rotation as well as the resulting 
Bell state. The bipartite entangled state |S')(|S'D) + \DS))/\/2 was produced 
with fidelity of 0.75. 



4-7 Process tomography 



Process tomography is a me t hod to characterize a quan tum mechanical evo- 



lution (IChuang and Nielsen! . 119971 : iPoyatos et al.l . 119971 ). Measurements are 



made to determine how an arbitrary input state, characterized by the density 
matrix pin, is transformed by the quantum process. The output density matrix 
Pout of the process can be expressed as 



2™-l 



Pout — Xij^iPin^j ■ 

i,j=0 



(28) 



Here Xij is the so-called process matrix, is the number of qubits and the 
operators Ai form a basis of the space of the 2^ x 2^ matrices. All rele- 
vant information on the quantum process is contained in the process matrix 
Xij- In the standard procedure, the 4^ separable states — {|0)i, |l)j, (|0)j + 
\l)i)/\/2, (|0)i + i\l)i)/\/2} for a single qubit — are prepared and then the 
output of the process is characterized each time with a full state tomography 
(3^ measurement settings). Inverting Eq. [28] yields the process matrix Xij- 

Such a process tomography has be en carried out for characterizing quantum 
gates in NMR ( Childs et al.l. 12001 ) and in linear-optics quantum computing 



( lO'Brien et al.l . 



20041 : iKiesel et al 



2OO5I ) . For ion traps a one-qubit process 
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(the teleportation of a aubit (IRiebe et all 120071 )) and two-qubit processes (a 
CNOT and its square (IRiebe et all l2006l )) have been characterized. 



Knowing the process matrix Xij fo^^ ^ill basic operations of a quantum com- 
puter is a very good basis for estimating the computer's performance. However, 
there are a few caveats: the number of necessary measurements to determine 
the process matrix Xij scales quite dramatically and thus it becomes quickly 
impractical to characterize processes with numerous qubits. Already for a 
four-qubit process, 20736 measuring settings would be required summing up 
to about 24 hours measurement time with the current parameters of the Inns- 
bruck experiment (100 experiments/setting, 25 repetitions/s). In addition, it 
is not clear that the subsequent application of two processes corresponds to 
the product of the process matrices. This assumption usually holds only if 
the relevant environment is time invariant, i.e. the process interacts only with 
a bath without memory. For example, for the Cirac-ZoUer gate, the phonon 
mode might keep the memory about the failure of an earlier gate operation 
and thus induce a failure of the next gate. Furthermore, the realization of the 
process can depend on whether it is executed at the beginning or the end of an 
algorithm. In particular, experiments are often triggered to the phase of the 
power line to reduce dephasing due to magnetic field fluctuations caused by 
50 Hz or 60 Hz noise and their multiples. Executing a certain gate operation a 
few milliseconds earlier or later within the experimental sequence leads easily 
to a change of the qubit resonance frequency of 100 Hz and the realization of 
the process becomes time dependent. 



Process tomography, as presented above, requires at least 4^3^ measurement 
settings and is thus quite costly. Above, we have restricted ourselves to the 
separable operators Aj. Using entangled auxiliary qubits and/or measuring 
in non-separable bases, the number of settings c ould be reduced, however , 
scales still exponentially in the qubit number N ( iMohseni and Lidan . 120061 ). 
Furthermore, most likely the total number of measurement runs has to be on 
the same order of magnitude as in the standard method to obtain a similar 
accuracy. Therefore, we conclude that a full quantum process tomography of 
a large quantum systems will be not practical. 



Finally, we note that there exist other approaches to estimate the fidelity of 
quantum processes. For instance, iKnill et al.l (120081 ) employ long sequences of 
randomly chosen gates. The main idea is that while the result of each gate 
sequence is known in an ideal implementation, noise leads to deviations from 
the expected results. Measuring the deviations, the average fidelity of the 
gate operations can be inferred. Choosing random sequences guarantees that 
the gate operation is investigated with various input states and in various 
combinations. 
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Fig. 31. General scheme of the Deutsch-Josza algorithm (from Guide et al. ( 20031 )). 
The upper line represents the qubit holding the input variable a which — if prepared 
in a logical eigenstate — does not change its value when Uf^ is called. The lower 
line holds the work qubit w. To this number the value of the function /„ is added 
modulo 2. The Hadamard rotations Ry^ (and Ry^) transfer the quantum bits into 
superpositions so that the inherent parallelism of quantum mechanics can be used. 

5 Algorithms with trapped ions 



5.1 Deutsch-Josza algorithm 



The Deutsch-Jo s za (DJ) algorithm de t ects t he parity of an unknown function 
(jPeutschl . Il989l : iNielsen and Chuang) . 12000). Concentrating on a single bit, 
there exist four different functions which map one (qu)bit with value a = 
{0, 1} onto another one. These functions can be divided into constant (even) 
(/i(a) = and /2(a) = 1) and balanced (odd) functions (/3(a) = a and 
/4(a) = NOT a). With a classical machine, it is necessary to call /„ at least 
twice to decide whether /„ is odd or even, i.e. one needs to calculate /„(0) 
and /n(l)- However, formulating the procedure quantum mechanically, the 
question whether is constant or balanced can be decided by calling it only 
once. 



In order to formulate the problem quantum mechanically, the functions /„ have 
to be generalized to take qubits as inputs. Within the framework of quantum 
mechanics all operations are unitary and therefore another qubit (the work 
or auxiliary qubit) is added to allow for non- invert ible functions /i and f2- 
Rephrasing the task, qubit \a) holds the input variable x while qubit \w) (the 
work qubit) will receive the result of the evaluation /„(a) plus the initial value 
w of qubit \w) to guarantee invertibility (see Fig. [31]). Thus, we define the 
unitary ?7/„ representing the implementation of the function acting on |ty)|a) 
with values w and a, respectively: 

UfJw)\a) = \Ua)(Bw)\a) . (29) 

Here, © denotes an addition modulo 2. 

The DJ-algorithm consists of the following steps (c.f. Fig. 1311) : 
(1) Initialize the system in the state |Oa)|lw). 
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Laser pulses 


fl=0 






/2 = 1 


Ry^. SWAP-1 
NOT„ 

SWAP Ry^, 


i?+(-^,0)i?+(^,(/.swap)i?+(;^,0) 

V V ^ V ^ 

ii^(§,0)ii^(^,§)ii^(f,7r) 

Vr)i?+(^, ^ + (^swap)^+(^, ^) 




Ry^ CNOT Ry^ 


i?+(^,0)fl+(7r,f)fl+(^,0)i?+(vr,f) 


fi = NOT X 


Ry^ 0-CNOT Ry^ 


i?^(7r,0) 
0)i?+(7r, f 0)i?+(7r, f ) 
ii^(7r,0) 



Table 1 

Laser pulses for the implementation of the algorithm inside the dashed box in 
Fig. [3T] {Rnj^^ U Ry^ ) on a single ion. (/^swap is given by arccos (cot^(7r/\/2)) . For the 
whole DJ-algorithm an i?'-^(7r/2, 0)-pulse just before and an i?*^(7r/2, 7r)-pulse after 
implementing Ry^Uf^Ry^ is applied. 



(2) Transfer the input qubit \a) into (|0) + |1)) /-\/2 and the work-qubit \w) 
into (|0) — |1)) /^/2 with Hadamard operations Ry. 

(3) Call the (unknown) function with these superimposed values by imple- 
menting Uf^. 

(4) Close the interferometer by applying an inverse Hadamard operation (Ry) 
on I a). 

(5) Read out the result in \a). 



The ion trap experiment used only a single ^°Ca'*'-ion ( j Guide et al.l . 120031 ). 
The internal state acted as the qubit \a) to hold the input variable for the 
function with the logical assignment |0) = \S), while the axial vibrational de- 
gree of freedom was used as the work qubit \w) (logical assignment |0) = |l)ax 
and |1) = |0)ax)- Thus, ground state cooling to 15", 0) initialized the system 
in |Oa)|lw), as required. A pecuharity of encoding a qubit within the ion's 
motional state is that one has to take care that the system does not leave 
the computational subspace {\S, 0), \D, 0), \S, 1), {\D, 1)}. In the experiment, 
this was achieved wi t h the composite pulse techniques described in Sec. 12.7.31 



fIChilds and Chuand . lioooh . Furthermore, single-qubit operations on the mo- 



tional degree of freedom had to be carried out. For this, the quantum informa- 
tion in the motional degree of freedom was swapped to the electronic degree 
of freedom such that ordinary carrier pulses could be used. Finally, the quan- 
tum information was swapped back to the motional degree of freedom. The 
Hadamard rotations Ry^ and Ry^, were absorbed into the definitions of the 
functions such that only for /2 this swapping of the quantum information had 
to be employed. 
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To decide on the class of the implemented function, only qubit \a) had to 
be measured. Finding \a) in |0) indicated that the function was even, finding 
|1) showed that the function was odd (see points at the end of the traces in 
Fig. [32]). For the functions /i, /s and the fidelity to identify the functions 
class with a single measurement exceeded 0.97, for /2 it was still above 0.9. 

In order to follow the evolution of |(l|a)|^, the pulse sequence was truncated 
at a certain time t and the qubit a was measured (see Fig. [32l) . Repeating this 
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Fig. 32. Traces for the implementation of the D J- algorithm (from I Guide et al 
The solid line is not a fit but a calculation based solely on the independently 
measured Rabi frequencies. The vertical dashed lines enfold the pulse sequences in 
Tab.H 

sequence for various times, it is possible to follow the algorithm through its 
evolution. This procedure helps to debug algorithms and makes sure that the 
desired algorithm is implemented. 



5.2 Teleportation 



In quantum teleportation, the state of a qubit is transferred from one physi- 
cal system to another one. This can be achieved with the following protocol 



( iBennett et al.l . 119931 ): first two parties share an entangled qubit pair. The 



quantum information contained in an additional qubit (the source qubit) can 
be transferred from one party to the other party by performing a Bell-state 
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Fig. 33. The teleportatio n algorithm' s quan tum circuit as implemented in the Inns- 
bruck experiment (from iRiebe et al.l (|2004l )). Double lines represent flow of classi- 
cal information, whereas single lines represent flow of quantum information. The 
gray lines indicate when a qubit is protected from detection light via so-called hid- 
ing-pulses. First, ions #2 and ^3 are entangled, creating the quantum link between 
the source region (ions #1 and #2) and the target ion (ion 7^3). Then, after some 
waiting time, the state to be teleported (on ion #1) is prepared via the unitary 
operation U^. A controlled phase gate together with detection via a photomultiplier 
tube (PMT) implements the Bell state measurement. 

measurement on the source qubit and one of the entangled qubits. To con- 
clude the transfer, the result of this Bell-state measurement is communicated 
via a classical channel to the receiver party and the receiver performs one of 
four rotations depending one the result of the Bell-state measurement. Thus 
it is possible to transfer the information content of a qubit by communicating 
two classical bits (the result of the Bell measurement) and using entangle- 
ment. Therefore, teleportation demonstrates a way to break down quantum 
information into a purely classical part and a quantum part. 



Another feature of teleportation is that it is not merely a simple transmission 
of a quantum state: it does not need a quantum channel to be open at the 
time the transfer is carried out. Instead it uses the non-local properties of 
quantum mechanics, established by a quantum channel prior to the generation 
of the state to be teleported. Once that link has been established, an unknown 
state can be transferred deterministically at any later time using classical 
communication only. Especially this feature w a s highhghted by the two ion 



trap teleportation experiments (iBarrett et al.l . I2004J : iRiebe et al.l . 1200J) by 



entangling the auxiliary and the target qubits before writing the quantum 
information into the source qubit. Thus these experiments demonstrate that 
unknown quantum information can be transferred on demand without using 
an active quantum channel. 

The teleportation circuit displave d in Fig. [33] is form ally equivalent to the 



one proposed by Bennett et al. (IBennett et al.l . Il993l ). but adapted to the 



Innsbruck ion -trap quantum pro cessor. Th e Innsbruck group reached fidelities 
of about 0.75 fiRiebe et al.l . l2004h a nd 0.83 dmebe et ahlj iooa), while the NIST 



group measured a fidelity of 0.78 (IBarrett et al.l . |2004| ). Teleportation based 



on a completely classical resource instead of a quantum entangled res ource 
yields a maximal possible fidelity of 0.667 (IMassar and Popescul . Il995l ). We 
note that to rule out out hidden variable theories, a fidelity in excess of 0.87 
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Fig. 34. Position of the ions within the segmente d trap during the exe cution of the 
NIST groups teleportation implementation (after Barrett et al. ( 20041 )). 



is required (iGisinl . Il996l ) . 



To emphasize the role of the shared entangled pair as a resource, in the Inns- 
bruck experiments, a delay between the creation of the Bell state and the state 
preparation of the input qubit was introduced. For waiting times of up to 20 ms 
(exceeding the time required for the teleportation by more than a factor of 
10) no significant decrease in the fidelity was observed. For longer waiting 
times, the measured heating of the ion crystal of less than 1 phonon/100 ms 
is expected to reduce the fidelity significantly, because the Cirac-Zoller gate 
requires the center-of-mass mode of the ion string to be in the ground state. 



The implementation of the NIST group (IBarrett et al.l . |200J) demonstrates 
how segmented traps facilit ate ion trap quant um computation. The authors 
use a segmented linear trap (IRowe et al.l . l2002l ). where the two DC-electrodes 
rails are each split into eight segments. Fig. [3l] shows the ion positions in 
each step during the teleportation procedure. By changing the potentials on 
the electrodes (top), the ion strings can be moved, split and merged. The 
teleportation algorithm was implemented in the following way (see Fig. 
first the leftmost and rightmost ions (auxiliary and target, respectively) were 
prepared in the Bell s tate — i|Tt)i.s u sing the geometric phase gate 

discussed in Sec. I^CT fiLeibfried et all . l2003bh . As the bus mode, the stretch 
mode was used. The center ion does not couple to the stretch mode and 
thus an effective two-qubit gate between the outer ions is implemented. For 
the experiments, it proved useful to transfer this Bell state into a singlet 
state (|\l/) = |ti)i3 — |iT))i3 as the singlet state remains invariant under global 
rotations allowing for an effective single-qubit addressing of the source ion #2. 
This was achieved by changing the relative position of the ions by varying the 
trap strength such that the subsequ e nt la ser pulses had the desired phases 
at the new ion positions (IRowe et al.l . l200ll ). To implement the teleportation, 
the Bell state measurement has to be carried out on the source ion #1 and 
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the auxiliary ion ^2. For this, the target ion #3 was first separated from the 
string in trap #6 and ions t^I and #2 were transferred back into trap ^5. Most 
importantly, the stretch mode of the two ions was still close to the ground state 
and the required rotation into the Bell basis could be directly implemented 
with the geometric phase gate (Sec. 12.6.41) without the need of a sympathetic 
cooling step (c.f. [6]). Then ion #1 and #2 were measured by transporting first 
the auxiliary ion into trap ^5, measuring there its fluorescence (see Sec. 12. 4p 
and pumping it into state ||). Next, ion #2 was also transported to trap #5 
and the total fluorescence was detected, revealing the state of the source ion. 
In a last step, ion #3 was transferred into trap #5 and conditioned on the 
result of the Bell measurement, the corresponding single-qubit operation was 
applied. 



5.3 Quantum error correction 



Classical computers use of error correction schemes intensively. It is to be ex- 
pected that quantum computers will employ error correction schemes as well. 
However, due to the continuous nature of quantum information, it might seem 
difficult to ap ply th e ideas of c lassical e rror c orrection to quantum comput- 
ers. However, IShoii ( 19951 ) and ISteand ( 19961 ) both found algorithms which 
correct errors by moving the errors from the quantum register to ancilla sys- 
tems. In these procedures, a logical qubit is encoded in a number of qubits 
such that the two-dimensional Hilbert space of a single qubit is embedded in 
a higher-dimensional space. Errors will then rotate the state vector out of the 
allowed subspace. Generalized measurements can project the system back to 
the allowed subspace and in case of small errors the original state is recovered . 



Quantum error codes we re implemented in NMR for the first time (I Cory et al. 



19981 : iKnill et all . l2001af ) 



Using trapped ions. IChiaverini et al.l (120041 ) implemented a rudimentary quan- 
tum error correction protocol. In these experiments, the authors encoded the 
(arbitrary) state of a source qubit in a superposition of two distinct three-qubit 
states (the primary qubit + two ancilla qubits), introduced controlled errors 
(spin flips only) on all three of them, before they decoded the state with the 
inverse operation used to encode the primary qubit. Small errors occurring on 
the encoded state rotated the state such that after decoding, the error could 
be corrected for. Read-out of the ancilla qubits provided the error syndrome, 
based on which the primary ion was returned in its original state. Using the 
language of quantum error correction, the stabilizer code \ZZX, ZXZ} (for 



1997 



an introduction to stabilz er codes see for instance Refs. (iGottesman 
Nielsen and Chuangi . |2000| )) was employed in these experiments. This encod- 
ing procedure was conveniently implemented by an entangling operation sim- 
ilarly to the ones discussed in Sec. 12.6.41 and used in Sec. 14.21 However, in this 
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particular instance, the three ions were placed in the standing wave such that 
they experienced the phases {(p,(p + 2/37r, + 4/37r} in the lattice. In this 
symmetric, configuration the ions felt no force if they were all in the same 
state. In all other logical basis states the total averaged force had the same 
absolute value and thus the same phase was acquired. Most notably, the heart 
of the algorithm (encoding and decoding) was executed only with global qubit 
operations, i.e. without individual addressing. Only for the preparation of the 
primary qubit, individual addressing was necessary and for detecting the error 
syndrome, selective state read-out was used. 

The NIST-group studied the performance of the error correction for the three 
different input states ^MTO\]) - iVOM\i) , V022\]) - i^/oJ8\i)} for ar- 
tificial error angles 9e applied to all qubits simultaneously while the qubit was 
protected by the encoding. Technical imperfections led to a fidelity of about 
0.8 even if no error was applied. Therefore the NIST group compared results 
where the error syndrome was used to results where the correction pulses were 
not used to correct the primary qubit. For the measurement eigenstate ||) and 
error angles 6e < 7r/2, the fidelity stayed close to 0.8, whereas in the uncor- 
rected case, the fidelity dropped to 0.5. For the two other superposition states, 
also a clear improvement over the uncorrected implementation was found. 

One of the biggest challenges in quantum information processing will be to 
improve the fidelity of an error correction algorithm such that it is below a 
fault-tolerant threshold. Furthermore, the qubit should never be left unpro- 
tected. This implies that the error correction has to be applied directly on 
the encoded qubit. Finally, it will be necessary to apply the error correction 
repeatedly and to extend the algorithm to correct for spin fiips as well as for 
phase fiips. 



5.4 Semiclassical quantum Fourier-transform 



The quantum Fourier transform is the final step in Shor's algorithm to factor 
large integers. It is designed to find the periodici t y of a quantum state (IShor , 



I994J: ICoppersmithl. 



I994J: lEkert and Jozsal . Il996l : (Nielsen and Chuand . l2000h 



Nielsen and Chuand (120001 ) show that the quantum Fourier transformation 
transforms an A^-qubit register in binary notation according to {ki G {0, 1}) 



k^kN-i ■ ■■k2ki) — > 
[(|0) + e2-[o-^il|l)) ® (|0) + e2-^[o.fefci]|^^) ^ . . . 

®(|0) +e2^^[°-'=^-i-'=2^'^]|l)) ® (|o) + e2^^[°-'=^^^-i-'=2'^i]|l))]/v^. (30) 
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Here [O.fc^fcjv-i ■ ■ ■ k2ki] stands for kN/2 + kN-i/4: + ■■■ + ^2/2^"^ + ki/2^ . 
Each qubit is rotated by i?*"('7r/2, — vr/2) and has acquired a phase shift condi- 
tional on other qubits. If the quantum Fourier-transform is the last step in a 
quantum algorithm, one can take advantag e of this structure and p erform the 
semiclassical quantum Fourier-transform ( Griffiths and Niu . 19961 ). Starting 
with the first qubit fci=i, first an Rf{n/2, —vr/2) pulse is applied to qubit fcj. 
Then qubit ki is measured and conditioned on the result exp{i7iaz /2^^~''^~^^^) 
rotations are carried out on qubits kj for all j > i. Note that a z rotation di- 
rectly preceding a measurement does not change the measurement result, such 
that the last z rotation before a measurement can be omitted. This procedure 
is then repeated for the next qubit i with i increased by one. 



Chiaverini et al.l (j2005bl ) implemented this procedure and tested it on a variety 
of separable and entangled three-qubit states. The implementation used a 
segmented trap, such that the state of one ion could be measured without 
destroying the quantum state of the others. In this way one after the other 
ion was measured and appropriate single-qubit operations were carried out 
conditioned on the measurement result. Four states (|001) + 1010) + ■ ■ --I- 1111), 
|001) + |011) + |101) + |111), |011) + |111), |111)) representing all the possible 
periods one, two, four and eight, respectively, of three bits were successfully 
tested. In addition, the entangled state aooi|001) + aoii|011) + aioo|100) -|- 
aiio|110) (I Oil = \aj\) with approximate period three was investigated. This 
state has the particular property that the result depends on the relative phases 
of the coefficients Oj which was confirmed in the experiments, too. 



5. 5 Entanglement purification 



High fidelity entanglement can be obtained from multiple entangled states of 
lower fidelity by a procedure called entanglement purification. In the context 
of fault-tolerant quantum computing, entanglement purificati on can alleviate 



thus t he stringent requirements for q uantum communication (IBennett et al. 



19961 : iGottesman and Chuangl . Il999l ): first many entangled states of a rela- 



tively moderate fidelity are created and shared between the two communicat- 
ing parties. Next, each party carries out local high fidelity quantum opera- 
tions and measurements before it communicates the measurement results to 
the other party (see Fig. [35]). Based on this, the parties can decide when the 
entanglement purification was successful and which states to keep for further 
use. This shared entanglement can be employed to perform a high fidelity state 
transfer between the two nodes via quantum teleportation (see Sec. 15. 2p . In 
this way, two quantum nodes can be linked with high fidelity even if the direct 
transfer of quantum information between them limits the fidelity of Bell pairs 
to just above 0.5. Together with quantum teleportation, this Bell state can be 
used for a high fidelity state transfer (see Sec. 15.21) . 
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Fig. 35. Schematics of a high-fidelity transfer of quantum information over a noisy 
quantum channel. First, a number of entangled states of moderate fidelity (> 0.5) is 
shared by both parties (quantum nodes). Next, local operations and measurements 
are used to create a single Bell state with high fidelity. Finally, this qubit pair is 
used to teleport quantum information from one node to the other. 

Reichle et al. demonstrated the entanglement purification procedure 

by purifying a two-atom Bell state out of two pairs of Bell states. First, the 
NIST-group used the geometric phase gate (see Sec. I2.6.4p to entangle four 
ions pairwise in a single step. To achieve this, the ions were placed such that 
ions #1 and #2 as well as ions #3 and #4 were spaced each by multiples of 
the standing wave (Ax = nX), while between ions #2 and #3 was a distance 
of (n ± 1/4)A (see Fig. 1361) . In this way the interaction between the two pairs 
vanished while within each pair the ions were entangled (ion #1 and #2 on one 
side and #3 and #4 on the other side). To implement the purifying procedure, 
a controlled-NOT operation has to be implemented on each quantum node, 
ion #1 and #3 on one side and ion #2 and #4 on the other side. For this, the 
axial confinement was changed such that there was a distance of nX between 
ions #1 and #3 as well as between ions #2 and #4, while the distance was 
(n± 1/4)A between ions belonging to different nodes. Thus another geometric 
phase gate executed two (local) phase-gate operations simultaneously. Mea- 
suring ion #1 and #2 (originally entangled with each other) in opposite states 
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■< ► 

X 

Fig. 36. Location of the ions inside the standing wave pattern for a given point in 
time during the phase gate operations used by the NIST group for entanglement 
purification. Panel a) shows the situation for entangling ion #1 and #2 on one side 
and for entangling ion ^^3 and #4 on the other side. The situation to carry out the 
two phase gates between ion #1 and ^^3 and between ion ^2 and #4 is illustrated 
in panel b). 

signaled success of the purification and the entanglement of the other entan- 
gled pair increased, while measuring them in the same state signaled failure 
of the protocol. In order to verify the entanglement between ions #3 and #4, 
the NIST group split these ions off from the other two ions prior to the mea- 
surement, such that the measurement did not destroy their coherence. In this 
way, a Bell state with a fidelity of 0.629 ± 0.0015 was obtained starting from 
two states of the form a|TT) + /3|ii) with a Bell-state fidelity of 0.614 ±0.0015. 



5.6 Quantum simulations 



Quantum simulations appear to be very appealing and especially fe asible ap- 
plicat ions of quantum computers. The idea was first put forward by iFeynman 
(Il982l ) and can be actually thought of as the birth of the field of quan tum 
information processing. Later Lloyd verified this suggestion (ILloydl . Il996l ) and 
together with Braunstein showed how a quantum system , with continuous vari- 
ables can be simulated efficiently (ILloyd and Braunsteinl . ll999l ). Quantum sim- 
ulations appear so attractive because they require only the interactions which 
are present in the system to be simulated. In addition, often only moderate 
manipulation fidelities are necessary to obtain meaningful results. 
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Leibfried et al.l (120021 ) show that the quantum dynamics of a spin-1/2 parti- 
cle in an arbitrary potential can be efficiently simulated with a single trapped 
ion. In the same publication the NIST-group simulates the action of an optical 
Mach-Zehnder interferometer and thus a first quantum simulation with a sin- 
gle ion. In this work, a beam splitter is implemented in the following way: the 
electronic state of the ion represents the number of photons in one of the in- 
coming mode^, while the motional state of the ion describes the state of one 
of the output modes. An i?+(7r/2,0) pulse on the blue sideband can transfer 
photons from one mode to the other and thus implements the beam splitter. 
Thus, a sequence of two i?"'^(7r/2)-pulses on the blue sideband, resembles the 
action of two beam splitters of a Mach-Zender interferometer. Furthermore, 
using second and third order sideband pulses (detuning A = nt^t, n = {2, 3}) a 
non-linear interferometer can be implemented where one photon in one mode 
can generate two and more photons in the other mode. This process is equiv- 
alent to parametric down conversi on — a quite impo r tant p rocess in quantum 



information science with photons (iTittel and Weihd . l200ll ) 



In the experiments, first a single ^Be"'"-ion was cooled to the motional ground 
state. Next, a 7r/2-sideband pulse of order n = 1, 2, 3 was applied on the axial 
motional sideband to create the state (|i,0) + |t,n + l))/-\/2. Then a change 
in trap frequency induced a phase shift between the two eigenstates. Finally, 
the phase shift was analyzed with another 7r/2-sideband pulse of order n. 
The phase shift acquired between the two n/2 pulses is proportional to the 
energy separation. Thus the fringes are more sensitive to the trap frequency 
change for larger n. This demonstrates the enhanced sensitivity of a non-linear 
interferometer as compared to a linear one. 

Quantum simulations with trapped io ns have also been inves tigated in the 



context of quantum phase transitions. iPorras and Cirad (120041 ) propose that 



each ion in the crystal represents a spin. Irradiating the string for example with 
off-resonant laser radiation produces in each ion a differential level shift which 
can be interpreted as a magnetic field acting on the spins. Interaction between 
the spins, i.e. an effective spin-spin interaction, can be simulated by coupling 
the ions to a collective vibrational mode. Thus a wide class of Hamiltonians 
can be investigated efficiently. In particular, it should be possible to observe 
quantum phase transitions in ion traps. 

A first step in this direction has already been taken by the Max Planck group 
in Garching by simulating the phase transition of a tiny quantum magnet 
consisting of two ^^Mg"*" ions from a parama gnetic to a ferromagnetic order 



( ISchaetz et al.l . l2004l : iFriedenauer et al.l . |2008| ). The spin-spin interaction was 



realized by the same laser-ion coupling used in geometric phase gates (see 
Sec. I2.6.4|) while the action of a magnetic field in the x direction was simulated 



Since only one photon is assumed at the input this description is complete. 
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by driving local spin-flips between the hyperfine states of the ions with an RF 
field of frequency cUqubit- The ground state of an ion with a magnetic field in x 
direction is given by |— ;>) = (||) + e*'^'i"''"*||))/-\/2- Tuning the ratio of the two 
simultaneously applied Hamiltonians adiabatically, the ground state |— of 
the parametric phase was transferred to the ground state |(TT+ii)/"\/2) of the 
ferromagnetic phase with a fidelity of 0.88. 



6 Shuttling and sympathetic cooling of ions 



Shuttling ions between various traps might rel ieve the requirements for scal- 



19981 : iKielpinski et al 



able ion trap quantum computing considerably (iWineland et al. 
2OO2I ). In accelerator experiments, shuttling of ions between different traps 



and re-cooling h as be en long established to slow down fast ions efficiently 



(iHerfurth et al.l . 120011 ) . Also single i ons h ave been transported reliably be- 



tween Penning traps ( Haffner et al. . 2000l ). However, for ion trap quantum 



computing, the requirements are more stringent: in particular, ion strings have 
to be separated, moved through junctions, recombined and the quantum in- 
formation must be preserved during all these operations. 



In this context, iRowe et al.l (|2002| ) demonstrated the reliable transport of sin- 
gle ^Be"'"-ions over 1.2 mm within tens of microseconds in a segmented ion 
trap. Moreover, they showed that the coherence of the hyperfine qubit was 
not affected by the transport. For this they first transfer the ion in a superpo- 
sition with a 7r/2-pulse, transported the ion and tested for the coherence. A 
contrast of 95.8% ±0.8% limited by magnetic field fluctuations was measured. 
In a second set of experiments, the NIST-group used a spin echo sequence to 
reduce the influence of magnetic field fluctuations. Here the ion (put again in 
a superposition of physical eigenstates with 7i/2 pulses) was transported back 
and forth. A spin-echo vr pulse was applied in the remote trap. The contrast 
of this measurement as well as of a control experiment without transport was 
almost 97%, limited by imperfections of the Raman pulses and off-resonant 
scattering from the P-level. These results demonstrated that the transport did 
not affect the coherence of the quantum state. 



Furthermore, iRowe et al.l (120021 ) investigated motional heating for various 
transport speeds. For this, the voltages on the trap electrodes were changed 
such that the trap minimum moved with time along a sin^ function. This 
ensures that not only the velocity but also the acceleration varies smoothly 
with time. They found that the ions can be transported between the two 
traps within 54 fis with hardly any observable motional heating (axial trap 
frequency 2.9 MHz, heating < 0.01 quanta/transport). Only for shorter trans- 
porting times heating of the ion motion was observed in agreement with a 
classical simulation. 
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In addition, iRowe et al.l (120021 ) demonstrate splitting of ion strings. However, 
their electrode structure was not optimized for this task as the size of the 
separation electrodes was too large. Therefore only success rates of 95% were 
achieved as well as relatively long separation times of about 10 ms were re- 
quired. Furthermore, the ions heated excessively by 140 ± 70 quanta due to 
the fact that during the separation the ion oscillation frequency was quite 
small. Later the NIST gro up demonstrated that these problems can be over- 
come (IBarrett et al.l . 120041 ) (c.f. Sec. 15.21) . Separation times of a few 100 /is as 
well as small motional heating have been achieved in these experiments. In 
particular, after splitting one ion from a three ion crystal, the stretch mode of 
the remaining two-ion crystal was still in the ground state, while the center- 
of-mass mode had acquired one motional quantum. A theoretical study of 
the splitting process and the required electrode structures can be found in 
Ref . iHome and Steand (120061 ) 



Another requirement for the proposal by Kielpinski et al. ( 2002ll is the trans- 



port t hrou gh junctions. First expe riments were carried out by iPearson et al 



(120061 ) and iHensinger et al.l (120061 ). In the former experiment, charged nano- 
particles were transported through a four-way crossing of a planar segmented 
trap. In the latter, Cd"*" ions were moved through a T-junction. Recently, the 
linear tr ansport of ions was studied within the framework of quantum me- 
chanics ( iReichle et al.. l2006ar). whi l e the non-adiabatic transport was investi- 



gated theoretica llv bvlSchulz et al.l (120061 ) and experimentally by lHuber et al. 



( 2008 ). Finally, Hucul et al.l ( 20081 ) analyzed the transport through various 
junction geometries (including T junctions) quantum mechanically. 



After the transport, the ion strings might have to be re-cooled such that 
the subsequent operations can be carried out with high fidelity. In order 
to achieve this while maintaining the coherence, a viable way seems to use 
the strong Coulomb coupling and to cool only a part of or even only one 
ion of an ion string. Such sympathe t ic cooling was demonstrated in yarious 



experiments (iDruUinger et al. 



1980; Larson et al.. 1986; Rohde et al.. 2001 



Hornekaer et al. . l200lh . Furthermore it has been established in various con 



texts ranging from precision measur ements (IRoth et al.l . |2005|) to electron cool- 
ing in accelerators (iMeshkovl . 119971 ) . 



Within the context of quantum information, iRohde et al.l (120011 ) cooled a two 
ion string to the motional ground state as required for many two-qubit propos- 
als. However, in these experiments two ions of the same species were used. As 
typically one of the qubit levels takes part in the cooling process, it is impor- 
tant to shield other nearby ions from the radiation to maintain the coherence 
of the qubits durin g the cooling proces s. This is very difficult if only one ion 
species is employed. iBlinov et al.l (120021 ) cool a ^^^Cd"*" ion via Doppler cooling 
of ^^^Cd"*". The isotope shift would help to preserve the coherence, however, it 
might be better to use a completely different ion species as done in the early 
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experiments by the NIST group. iBarrett et al.l (120031 ) cooled Be"*" with Mg"*" 
to the motional ground state (vice versa only to the Doppler limit) and in 
separate measurements showed that even on time scales of 30 ms a qubit en- 
coded in the hyperfine manifold of ^Be^ is n ot affected by the strong cooling 
light for Mg"*". Finally, ISchmidt et al.l (120051 ) took the idea one step further 
and cooled and detected the internal state of A1+ via Be"*". 



While all prerequisites for quantum computing by shuttling ion strings have 
now been demonstrated in separate experiments, the combination of shuttling 
ions, splitting and re-cooling the ion strings in the same experiment and at the 
same time preserving the quantum information has yet to be accomplished. 



7 New trap developments 



Parallel to the efforts to shuttle and split ion strings, in particular the NIST 
group has put quite some effort into developing new traps manufactured by 
microfabrication techniques to build a medium sized quantum computer. Mi- 
crofabrication techniques allow for complicated and precise electrode struc- 
tures. Trap sizes (measured as the distance between two RF electrodes) range 
from about 200 fim down to a few tens of yum. While in the beginning mainly 
three-dimensional designs were b uilt using two sub s trates, recently a new, pla- 
nar trap design was invented by IChiaverini et al.l (j2005al ) . In such traps the 
electrodes are aranged all in a plane such that a single substrate suffices to 
mount the electrodes. Typically there are five (four) electrode groups in such 
a linear surface trap (see Fig. [371) : on the center electrode, ground or a small 
DC potential is applied. The two neighboring electrodes receive an RF po- 
tential which provides the radial confinement. The outer electrodes are often 
structured and various DC potentials provide the axial confinement. The ions 
are trapped in the center of the field lines in Fig. l37b where the RF field 
vanishes. The asymmetry of the RF electrodes tilts the principle axes of the 
quadrupole field such that a laser beam travehng parallel to the trap surface 
has a projection along both radial motional degrees of freedom and all modes 
can be cooled satisfactory. Two trap designs were suc cessfully used by the 



NIST group, one based on gold on a qu artz substrate (jSeidelin et al.l . l2006l ) 



and the other one based solely on silicon (iBritton et aD . l2006h . Ions have been 
trapped as close as 40 fim to the sur f ace w i th still reasoii a ble h eating rates 



of a few phonons/ms (ISeidelin et al.l . l2006l : lEpstein et al.l . l2007l ). The MIT 



group developed a tr ap based on printed circuit board fabrication techniques 
( Pearson et al.l. 120061). The v also built traps using silver on a quartz substrate 
(ILabaziewicz et al.l . l2008al ). For trap sizes between 75 /im and 150 /im (ion- 
substrate distance), heating rates between 2 and 20 phonons/s were measured 
at 4 Kelvin (see Sec. 
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Fig. 37. Electrode configuration of a planar trap. Fig. [37b shows the a three-dimen- 
sional view of the asymmetric electrode arrangement, while Fig. 137b shows a side 
view with some of the electric field lines. 

8 Future challenges and prospects for ion trap quantum computing 



In order to achieve universal quantum computing, the algorithms have to 
be implemented in a fault-tolerant way. It is commonly accepted that this 
requires quantum error correction. Therefore, currently one of the most im- 
portant goals is to implement quantum error correction repeatedly with high 
fidelity to prolong coherence times and to correct for errors induced by the 
gate operations. The largest obstacle to perform a successful quantum error 
correction protocol seems to be the limited fidelity of the operations. The cur- 
rent state of the art for the control in ion trap quantum computing can be 
summarized as follows: 



The qubit coherence times are one or two orders of magnitude longer than 
the basic (gate) operations. In specific cases coherence times longer by more 
than five orders of magnitude the gate time are available (see Sec. I3.1.2p . 
In most current experiments, motional decoherence is not a problem. In 
addition, it can be further suppressed with cooling of the trap electrodes 
(see Sec. 13.21) . 

Initialization accuracies are on the order of 0.999 as discussed in Sec. 12.41 
Most likely they can be improved further if necessary. 
Single qubit oper ation can be carried out with fidelities exceeding 0.995 
(iKnill et al.U2008l ). If required, further improvements are possible with more 
stable laser fields at the ion positions. 

Implementations of two-qubit gate operations achieve fidelities of about 0.9- 
0.99. Depending on the gate type, various sources limit the fidelity. Errors 
are caused by off-resonant scattering, imperfect addressing of individual 
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Fig. 38. Prog r ess i n reducing the error rate of two-qubit gates (taken from 
Ref. iBenhelml (j2008l )) and in increasing the number of entangled ions. Open cir- 
cles represent experiments using two-qubit gate operations with global addressing 
while diamonds show results based on individual addressing. The performance is 
measured in terms of the infidelity of produced Bell states. The stars mark the 
largest number of entangled ions obtained at that time. Numbers below the ref- 
erence indicate the number of trap cycles required for the operation. Dashed and 
dotted lines indicate the trends. 



qubits, insufficient cooling, laser frequency and intensity noise. 

The read-out of a single qubit can be performed with a fidelity of 0.999. 

Further improvements seem possible (see Sec. 12.41) . 

Ion strings can be shuttled, split and merged (see Sec. E]) with high fidelity 
and small decoherence. 



Currently, two-qubit gate operations seem to be the main limiting factor and 
receive therefore most attention both from experimenters and theoreticians. 
Figure [38] shows the pro g ress of the fidelity achieved in the last decade. Most 
notably IBenhelm et al.l (j2008bl ) demonstrate two-qubit gate fidelities high 
enough to allow in princi ple fau l t tole rant quantum computation according 
to the scheme proposed by iKnilj (120051 ). In addition, one must always keep in 
mind that all these requirements have to be met at the same time. Further- 
more, so me emphasis should be given to parallel processing quantum infor- 
mation (jSteand . |2004| ). However, also initialization of the necessary ancillas, 
read-out, coherence times and the particular layout and the attainable degree 



81 



of parallelization are important. 



Both, analytical and numerical results, indicate that operational fidelities on 
the order of 0.9999/operation seem to be sufficient to achieve fault tolerance 
(the so-ca lled fau l t-tole rant threshold), provided certain other criteria can be 
met, too (Steane, 20041 ): specific errors, error propagation, the allowed over- 
head, specific requirements, the amount of possible parallelization, amongst 
others have to be considered to get a full grasp on the situation at hand. Thus 
the co ncept of thresholds is oversimplifying the situation. For instance, iKnill 
( 120051 ) published numerical results which indicate that even error rates on the 
order of 10~^ are permitted, however with a huge overhead of 10^ physical 
qubits for one logical qubit. It seems reasonable that every operation in a 
quantum computer should be treated with the same meticulous attention and 
be implemented as perfect as possible to achieve fault tolerance while keeping 
the overhead as small as possible. 



Interestingly, state transfer between interconnected ion trap quantum com- 
puters at the fault-toler a nt level is not necessary a s discussed in Sec. 15.51 
( Gottesman and Chuang . 1999 : Reichle et al. . 2006bl ): quantum information 
can be teleported deterministically between two locations using a purified en- 
tangled Bell-state. Entanglement between distant ions can, for example, be 
generated by splitting an entangled two-ion string and transporting one of 
the ions or by interfering fiu orescence light frorn two ions on a beam sp litter 
as demonstrated a l ready by iMaunz et al.l (120071 ) , iMoehring et al.l (120071 ) and 
Matsukevich et al.l (l2008h . 



Another important issue is the speed of the operations. Having in mind uni- 
versal quantum computing to outperform classical computers, e.g. in factor- 
ing large integers, billions of operations have to be carried out. Thus, the 
current typical time scales for the basic operations of a few hundreds of /xs 
seem simply too slow for factoring large numbers even if the operations are 
carried out to a large extent in parallel. There exis t proposals for gate opera- 



tions which are faster t han the trapping frequencies (jGarcia-RipoU et al.l . 12003 



20051 : IZhu et al.l . 120061 ) . however, there are other bottlenecks such as read-out 



and ion string separation which might slow down the processor speed. 



From this perspective, it seems attractive to work on hybrid devices where e.g. 
quantum information is stored in trapped ions (incl. error corre ction) and most 
of the processing is i mplemented e.g. with Josephson junctions (iMakhlin et al. 



oi tne processmg is i mplemented e.g. witn Josepnson j unctions (llviaKnun et al. , 
200ll : ISteffen et all l2006l : IPlantenberg et all . l2007h which operate at speeds 
roughly three orders of magnitude faster than current ion trapping approaches. 
Before this can happen, however, tremendous difficulties have to be solved in 
transferring quantum information between the two systems within the short 
coherence times of the Josephson junction qubits. Currently, t here exist a 



few proposals to couple ion-trap and Josephson junction qubits (iTian et al. 
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2004 l2005l : IS0rensen et al.l . |2004| ). but so far almost no experimental results 



were achievable. Furthermore, within the field of Josephson junction quantum 
information processing, there are some open challenges before advantage can 
be taken of the speed of this system. The biggest of these seems to be that 
coherence times of only a few /xs have been achieved so far. 



The original proposal of Cirac and ZoUer (jCirac and ZoUerl . Il995l ) is scalable 
in the sense that it does not require exponentially many resources with an 
increasing number of qubits. Indeed, Fig. [3H] shows that in the last decade 
larger and larger ion strings have been entangled. However, it seems imprac- 
tical to construct a device for a large number of qubits by trapping all ions in 
the same trap because it gets more and more difficult to obtain the required 
strong radial confinements to work with linear ion strings at reasonably high 
axial trapping frequencies. Furthermore, the mode structure of the ion crystal 
gets more complicated with more ions as well as the speed of the sideband 
operations is reduced with the larger mass of the crystal. Up to date there 
are a couple of routes known which potentially ease these technological chal- 
lenges. Almost all of them are based on distributi ng the ions across different 
traps and to intercon nect these traps via photons (Cirac et al.l, 
conducting strip li nes (ITian et al.l.l2004l:lHeinzen and Winelandl . 
via auxiliary ions (jCirac and ZoUen . 12000) • 



1997h . super- 
or even 



The currently most advanced proce dure, however, i s to merge and shuttle 



small ion strings in segmented traps (IWineland et al.l . Il998l : iKielpinski et al. 



2OO2I ). Those ideas have been studied in detail and seem to offer a practicable 



and viable way to scale ion trap quantum computers. Currently, big efforts are 
under way in realizing this architecture (c.f. Sec. [H]). Major challenges are the 
fabrication of such complex small ion traps combining high flexibility of ion 
movement (junctions), low motional heating rates and high trap frequencies. 
It is also strongly desirable to integrate the control electronics and optics on 
such ion trap devices. 



In summary, the basic requirements for a general purpose quantum computing 
device with trapped ions have been demonstrated and no fundamental road 
block is in sight. However, building such a device is extremely challenging. 
Especially, the stringent requirements for fault tolerance and for scalability 
to many thousands of qubits pose huge difficulties. However, reaching a good 
control over a reasonable number of qubits seems feasible in the next decade 
and might be of quite some interest: already with about forty qubits, phys- 
ical systems can be simulated which are intractable with current computing 
technology. It remains to be seen whether and how the dream of universal 
quantum computing can be implemented with trapped ions. 
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